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LANCHAD 


REMOTE SENSING OF THE N’DJAMENA ZONE AND 
THE LO GONE-CHARI CONFLUENCE 

Diachronic and Multispectral Cartography at 
1/50 3 000 and 1/100,000 

Claude Bardinet*, Jean-Marie Monget** 


- • FOREWORD 

For an understanding of the Earth and for an inventory of its 
resources , remote sensing- from, satellites offers a new means, the 
possibilities of which one is only -'beginning to perceive, and which 
supplement rather than substitute for the previous modes of inventory. 
Satellites are furnishing, in increasing torrents without cease, data 
on the totality of our planet. For the Republic of Chad, LANDSAT 1,' 
LANDSAT 2, the NOAA series, METEOSAT, and doubtless soon SPOT, have 
furnished, are furnishing, and will furnish a mass of multispectral 
data-thermographics, the exploitation of which has today widely begun 
at least under the heading experimental; but the time is near in 
which such exploitation will comprise a necessary foundation for all 
undertakings for economic development and management of the area. 

The LANDSAT data in particular appear well suited to the geo- 
graphic study of the Republic of Chad because they allow the estab- 
lishment of cartography at the 1 / 50,000 scale more detailed than that 
which has been available until now. Their spatial resolution also 
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permits a good recognition of the regional boundaries which 
delineate homogeneous characteristics over fairly large surface 
areas. Finally, by their a grasp of 

the variations in the vegetative cover and. of the water areas. 

In addition to these general advantages, remote sensing offers 
another — that of the rapidity of data acquisition. To this must 
be added that of .the processing and distribution of the results, 
because states prefer the preparation of true living atlases, con- 
tinuously updated, ‘ basic to the undertaking of economic development, 
rather than an historic atlas based on conditions already past in 
which a retrospective understanding is much less useful. 

This section presents the first results of the LANCHAD remote 
sensing program. One will see that this new tool for geographic 
research bases its interpretations on a necessary association of 
the knowledge of the terrain and the processing of the data by 
computer. 


Fernand Verger 

Professor, Eeole Nprmale Superieure 
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INTRODUCTION 


The LANCHAD Atlas is a part of the regional atlas of Chad, the /3 
work on which began in 1977 at the request of the Ministry for 
Economics, Planning, and Transportation of Chad made to the Fund for 
Aid and Cooperation. 

These regional atlases have been conceived principally to be a 
collective work by the following entities: Planning and Development 

Authority of Chad, the University of Chad, O.R.S.T.O.M. , and the 
University of Paris I. 


The scientific committee for this atlas was set up by Jean 
Cabot, Professor at the University of Paris VIII, Daniel Noin, Pro- 
fessor at the University of Paris I, and Gilles Sautter, Professor 
at the University of Paris I. 


The section "Atlas LANCHAD" has been directed by Fernand Verger, 
Professor at the Montrouge Teachers College. The present publica- 
tion LANCHAD has been carried out by Jean-Marie Monget, Director of 
the Center for Remote Sensing and Analysis of the Natural Environ- 
ment at the National College of Mines in Paris and by Claude Bar- 
dinet. Instructor at the Laboratory for Cartography at the Univer- 
sity of Chad. On September 30 3 1977 3 the team which had collabo- 
rated on the LANCHAD project also included Michel Albuisson and 
Michel Poisson of the Center for Remote Sensing and Analysis of the 
Natural Environment. 

The automated cartography was done at C.I.R.C.E. of C.N.R.S. 
in collaboration with Gerard Joly. 
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I Anne Le Fur (University of Paris I - E.S.C.G.) put together the 

I dummy press layouts for the press at C.E.R.C.G. of C.N.R.S. J 

}■ i 

f i 

I i 

j The financing of this first phase of LANCHAD was provided by the ; 

I F.A.C. within the framework of an inter university accord between the \ 

| University of Chad and the University of Paris I. j 

j ■ 

\ The methodology used is that perfected in the framework of the I 

| Fralit team (RCP 353 of the C.N.R.S.) by the C.T.A.M.N. The thematic j 

S cartography of the regions is here derived by a numerical treatment 

| 

i of the multispectral data obtained using satellite remote sensing. ; 

l The relationships existing between the nature of the objects on the > 

t ground and their multispectral signature permit a computer-assisted 

exploration and cartography of the regions. In the spectral regime 
| of the LANDS AT satellite (visible and near infrared) the dominant 

I problems are linked to chlorophyll, to water and to the nature of the 

| soil, while for the NOAA satellites (thermal IR), the phenomena ob- 

I served are more physical: geology, geomorphology, climatology. 
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CHAPTER 1 - HISTORY. OP THE LANCHAD ATLAS 


The LANCHAD atlas project was carried out within the framework 
of the interuniversity accord associating the University of Chad and 
the University of Paris I in the development of, teaching of, and 
research in geographical cartography. 

Geographical cartography being, in this choice, principally 
considered to allow the creation of maps useful for the management 
of the area, the thought has led to a regional type of atlas which 
has been lacking for the planner in our times. The coincidence of 
this objective of university research and of the request expressed by 
the Chad planners who wanted to quickly arrange a cartographic in- 
ventory of Chad in the human, economic, and physical regimes, led to 
the organization of a symposium* at N’Djamena on the subject of 
’’Regional Atlases of Chad and Regionalization of the Plan”. Of the 
different subjects covered in the course of this symposium, we will . 
here recall the one entitled: "The LANCHAD Project: Remote Sensing 

and Cartography of Chad”. 

The term LANCHAD, being indeed the contraction of ’’LANDSAT 
Imagery of Chad" for "Land Use in Chad", is the term chosen to empha- 
size that one intends to use the LANDSAT recordings for practical 
studies on the management of Chad territory. 

The first work in actual application of remote sensing tech- 
niques to the imagery of Chad was introduced in September 1976, 
after a one-year study period in the course of which different 
studies on the handling of the imagery were successively considered 


This symposium was convened at the University of Chad on January 
14-15, 1977, consisting of the geographers from the University of 
Paris, J. Cabot, D. Noin, B. Rouleau, and G. Sautter, the university 
staffs, planners of Chad and of researchers; the work will be pub- 
lished by the University of Paris I. C. Bardinet read the paper 
"Project LANCHAD: Dichronic Analysis of the LANDSAT Remote Sensing 

Imagery of Chad in the N’Djamena Zone at 1/50,000 and 1/J00,000. 


by C. Bardinet, in taking account of the experience acquired at the 
School of Advanced Applied Studies and the School of Advanced Studies 
in Social Science,, and by the teams ib^dii^it together, one around 
F. Verger (Fralit) and the other around J.-P. Gilg and G. Sautter. 

It is finally the methodology of the numerical processing pro- 
gram for tapes developed in the framework of the Fralit effort* which 
was employed to^§|^j|;:pip;.maps presented in this first LANCHAD pub- 
lication. The work method is called that of spectral association /5 

and taxonomical association after analysis. J.-M. Monget led the 
unmonitored analysis 3 and C. Bardinet — the monitored. The FRACAM** 
system used was developed at the C.T.A.M.N. of the National College 
of Mines in Paris. 

The test zone presented in these first sections is that of 
IvPDjamena (formerly Fort Lamy), capital of Chad, at the regional 
scales of 1 / 50 , 000 , 1/100,000, 1/200,000 and 1/500,000. 


* 

The Fralit Team - "Remote Sensing of the Oceanic Littoral of France n . 
’ /Cbiiectlhh de/iL f Ecole Normals Superieure de Jeunes Filles , No. 11, 
1977 3 310 pp. 

J.-M. Monget - .Automated Classification of Multispectral Data 
Using Harmonic Analysis. The CLAMS System. Bulletin of the French 
Society of Photogramme try , 1967 - 62, pp. 15-23. 


6 


f 


_ 1 

.-j 

CHAPTER 2 - THE GEOGRAPHICAL SETTING | 

The N ! Djamena zone (12° 6 north and 15° east for the center of /6 • 

the cluster) experiences a climate of the Sahel type with two sea- ] 

sons — a wet season from June to September and a dry season from 
October to May. Rain measurements are on the order to 600 mm per 
year with strong variations during the years ( 602.8 mm in 1972 , 

314.7 mm in 1973 3 and 990.1 mm in 1959 3 according to Tobias ] 

ORSTOM ;, 1976) . I 

1 

The atmospheric humidity varies from 80% in August to 21% in 

i 

March. Northeast winds predominate. The mean annual temperature 

is 28°, with a minimum of 14° in January and a maximum of 4l° in 

April (27-year averages established by ORSTOM). -Diurnal variations 

are important in all seasons, and can reach 20° in January and 8° J 

in August. The annual evaporative transpiration was estimated by 

Tobias (1976) at 1900 mm. ? 

2.1. The soils i 

A study of the soils (Pias and Barberry, ORSTOM, 1964) is pos- 
sible from the 1/200,000 map. 

2.1.1. The right bank of the Chari-N^jarema zone 

One observes the following dominances: 

Halomorphic alkaline soils of generally temporary hydromorphy: 

• cultivated soils in topographical depressions, related to 
argilaceous materials (black tropical argile); 

° soils of clay and pseudo-clay in patches and concretions, | 

related to recent sand-silt and clay-silt alluviations J 

(lateral ridging of the wadis and the drainages). f 


2.1.2. The left bank of the Chari and the 
Logone -Chari confluence 


One observes the following dominances : 

• cultivated soils in topographical depressions, soils of clay 
related to argilaceous material, here and there in chalky 
nodules; 

• hydromorphic mineral soils generally of semi -permanent hydro- 
mo rphy : 

- clay and pseudo-clay soils in patches and concretions, 
related to recent sand-silt and clay-silt alluviations 
(depressions and terracing of the Chari and the Logone), 
generally of semipermanent hydromorphology; 

- hydromorphy general, or of temporary depth, clay and 
pseudo-clay in patches and concretions, related to recent 
sand-silt or clay-silt, lateral ridging in wadis and 
drainage ; 

• complex hydromorphic-deposition (interfluvial) soils with /7 

alternate : 

- raw hydromorphic soils of sand or sandy argiles; 

- soils of argilaceous material here and there in 
chalky nodules. 

One observes a very great heterogeneity of materials at the 
outcroppings . 

The floodable areas are either of the recent argilaceous 
lacustrian series or of the recent sandy series of fluvial origin. 

The features above the water are either rims of the bank or of 
water courses or the contours of the old sandy series. The pumping 
of underground water slightly entrains the presence of salts. 


2.1.3. The meanders in the Chari and the Lo gone- 
in the N : Djamena zone 

They are characterized by: 

- the recent sandy series of features of quartzite sand, 

• argilaeeous stratification 

• the old argilaeeous -sandy series of chalky nodules (banks of 
the Chari), 

• the older alluvial series with current fine textured silt, 
silt-clays, and clay-silts (ridging of the banks of the 
rivers and tributaries, terraces and principal layers). 

The nature of the material on the surface determines in large 
part the level of reflectance recorded using remote sensing. Tobias 
(1976) determines in this zone four principle types of soils, either 

" hydromorphic soils on the rims of the bank, silty, of an 
excess of water (useable reserve drained off dry in 24 days), 
fine textured, compact, slight permeability; 

• ferruginous tropical soils of deep sandy material, fine 
textured, sandy on the surface, argilaeeous at depth, compact 
and not much permeability (useable water reserve drains very 
dry in 4 days); 

• cultivated soils "Berbere (sorghum) and ’’Karal", clay content 
elevated at the surface, heavy soils, fine constituents, 
argilaeeous minerals swelling with the formation of mud when 
water soaked and cracking with dessication in the dry 'state, 
slight permeability (useable reserve of water drains off dry 
in 10 days) ; 

• halomorphic soils or nagas, horizon sandy on the surface, of 
variable thickness, impermeable at the argilaeeous level, 
retention of water--. slight (reserve of useable water drains 
off in 5 days in dry season). 


2.2. The vegetation 
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In this zone, Pias ( 1968 ) observes six type of vegetation. 

2.2.1. Sandy wooded savanna (north of N 1 D,j amena ) s part ially 

on the rims of old or current watercourses and partially hydromorphic 
raw alkaline soils, dominated by Acacia senegalensis , Balanites 
aegyptiaca. Acacia seyal, second growth of Hyphaena thebaica. 

2.2.2. Vegetation on bottomlands and on nagas soils 

The nagas soils are often dispersed and surrounded by flooded 
lowlands. The lowlands are wooded as a function of soil texture. 

The vegetation is composed of: 

• small islands of dense vegetation, tall trees, Zizyphus 
mauretania, Acacia seyal, Anogeissus leiocarpus and gramineous 

T 7 . 1 v’lAarpe t s 'jo f tall b ear d gras s ; 

• open savanna on black tropical argils. Acacia seyal, and tall 
beard grass:; 

• Acacia seyal and Acacia scorpioides on the ponds marking out 
. . the nagas soils. 

2.2.3. The undergrowth 


Borassus aethiopum, Ficus, Faidherbia labida. Balanites aegyp- 
tiaca, Calotropis procera, on light textured soils where the water 
table is near the surface. 

2.2.4. The niarshy plain and the flood plains 
of the Logone 


E 

-r 


This type includes the flooded areas of the rivers during the 

rainy season, the interior basins, the edges of the wadis, the per- 
manent marshes. Here, flooding is variable from August to December 
The vegetation includes: 
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• in general, a gramineous carpet of beard grasses on hydro- 
morphic argilaceous soils ; 

• in the low lying areas forming puddles in the dry season, 
a carpet of Nymphea lotus; 

• on the edge of the flood plains or on knolls 5-6 meters 
above the water, Fisuc, Borassus aethiopum, Faidherbia 
albida, Hyphaene thebaic a, Tamar indus indica, Kigelia 
africana. . . 

2.2.5. The forest arcade, densely (‘wooded savanna : 

This is localized on the rim of the Chari and in part on that 
of the Lpgone 7 . The alluvial soils of recent origin, with an acid 
pH, of argilaceous silt texture support: 

• large trees, all types of acadia, Mitragyna Africana, Hyph 
Hyphaene Thebaic a. Ficus Livers, Zizyphus jptjuba; 

■ thickets of shrubs and thorns, Zizyphus Mauretania, acacia, 

C ombre turn ... 

• ih the low lying areas a pure population of Acacia seyal 
(including in the courses of the argilaceous wadis), 
along the banks, spiny bushes, mimosa asperata... 

2.2.6. The open savanna of black tropical argile 

This is a savanna of little or of limited flooding. One finds 
Acacia Seyal and Mitragyna Africana. 

2.3. Cultivation 

In the period of observation, at the beginning of the dry sea- 
son from October to February, there is principally seen: 


• cultivation of the sorghum called n berbere u in the low lands ; 

• truck farming of the more recent alluvial series of fine 
silty loam with a little mica content (Tobias), principally 
on the rims of the river banks and on the terraces in pro- 
portion to the water table. 



CHAPTER 3. THE INFORMATION METHODS EMPLOYED 


They comprise three information system of different character- /10 
istics , each belonging to a particular administrative entity, namely: 

- the Information Processing and Industrial Statistical Service 
(STISI) of the Ministry of Industry, Commerce, and Handi- 
crafts (MICA) equipped with a CII-IRIS 80 digital processing 
system; 

- the National College of Mines of Paris (ENSMP) equipped with 
an HP-2100 system; 

- the Interregional Center for Computers (CIRCE) of the National 
Center for Scientific Research, equipped with an IBM 370-175 
digital processing system. 

Each of these information systems -was used in a particular mode 
to support the programs developed either by CTAMN or by the Labora- 
tory for Geography of the ENSJF (F. Verger). 


3.1. The CII-IRIS 80 system of STISI 


The Computer Center of STISI is equipped with a CII-IRIS , 

80 biprocessor with a one-million 32 -bit word memory functioning 
with a variety of peripherals, such as bulk storage disks (200 million 
8-bit bytes) and 9-track, 1600 BPI tape decks. 


The ENSMP is connected to this computer by virtue of the parti- 
cipation of its parent ministry (MICA) which lacks information sys- 
tems . Thus, the CTAMN has the possibility of using this system in 
one of its two modes: operation by blocks share or by time sharing. 


For the remote block sharing , the ENSMP at Sophia Antipolis is 
equipped with a larger terminal SEMS-MITRA (Figure 3-A) connected to 
the STISI computing center at Voluceau by a line of 9600 baud capa- 
city. The large terminal is composed of the following elements: 
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- a SEMS-MITRA minicomputer 15 with a memory of 16K 16-bit 
words ; 

~ a card readout; 

- a printer; 

- a command console. 

This mode of operation is used by the DUERB 3 LICLA 3 AND DUCLA 
programs developed at CTAMN : 

- the DUERB program transfers the data received from the EROS 
DATA CENTER on magnetic tape to bulk storage disk, thereby 
permitting a subsequent use of the data in time sharing. 

- The LICLA program creates a printed version (listing) of the /12 
classification results using a variable key. 

- The DUCLA program enables a transfer of the classification 
results to magnetic tape thus permitting them to be used on 
the IBM-370 system at CIRCE. 

For shared time operation , the CTAMN at Sophia Antipolis is 
equipped with a small TEKTRONIX terminal (Figure 3-A) connected to 
the computing center of STISI at Voluceau by a line of 200 bauds 
capacity. This small terminal is composed of the following elements: 

- a TEKTRONIX 4006 display console 3 

- a hard-copy system on TEKTRONIX thermosensitive paper. 

This mode of -interactive operation is used by the FRACAM system 
of classification developed at the CTAMN. We will review in more 
detail (in Chapter 5) the different programs which comprise this 
system. They simply specify that the classification records are 
able to be produced at a scale that varies continuously from l/ 25 3 000 
to 1/1 5 000 3 000. 

3-2. The HP-2100 system at CTAMN 

For their own needs 3 the Center for Remote Sensing and Analysis 
of the Natural Environment (CTAMN) is equipped with an autonomous 
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modes of operation of the 








processing system constructed around the Hewlett-Packard 2100 mini- 
computer (Plate A). It comprises the following elements: 

- a minicomputer HP-2100 with a 16K l6-bit word memory; 

- two storage disks, each of 10 M 8-bit byte capacity; 

- a tape deck of 9 tracks-800 BPI; 

- a TEKTRONIX 4013 display console; 

- a VERSATEC 1200A printer-plotter. 

This processing system is used by the DUPSTR and VISTR programs 
developed at CTAMN. 

- The DUPSTR program permits the transfer to disk storage of 
the data received from the EROS DATA CENTER on magnetic tape, 
this choice perhaps associated with the spectral criteria 
utilizing one of the channels (MSS 7 for example). 

- The VISTR program creates a printout of selected data in the 

form of a black and white document of ten gray scale levels 
using the VERSATEC printer-plotter. These printouts can be 
created at different scales, respectively: 1/1,000,000, 

v 1/500,000^ 1/400, OOti, , and 1/100,000. 

The base. Used is elthefc white paper or transparent tracing 
; paper. 


This method of cartography has been used within the framework /l4 
of the LANCHAD project to print out the data from the LANDSAT and 
NO A A 5 satellites. 

3 ■ 3 . The IBM 370 system at CIRCE 

The CIRCE Computer Center is equipped with an IBM 370-175 
system of great capacity, functioning with the classical peripherals 
(9T-1600 BPI tapes and disks), as well as special peripherals among 
which is an incremental plotter of the BENSON type. 
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System of remote sensing data processing 
using a minicomputer (processing center 
of the CTAMN, College of Mines) 

Plate A. Different elements of the 
automated processing system used 
for the LANCHAD project 

(Plate continued on following page) 

This apparatus permits the creation of color maps from the 
classification output using the FRACARTO cartographic system devel- 
oped in the Laboratory for Geography of the ENSJF. 

These maps can be generated at all scales running from 
1/1,000,000 to 1/50,000, either directly as polychrome maps on 
opaque paper, or as a series of prints corresponding to each color 
of the printed map, each of such prints being plotted in black ink 
on transparent paper of the first quality (mylar). 



REPRODUCIBILITY OF THE 

ORIGINAL PAGE IS POOR 


Plate A. (continued) 







3-4. An example of processing (PGR 3) 

In order to familiarize the reader with the times of operation 
pertaining to remote sensing data processing, we present an example 
of this in the case of data processing for the N'Djamena area on 
January 12, 1973 (FOR 3). 

3.4.1. Utilization of the CII-IRIS 80 system 

- selection of a test zone in the neighborhood of 40,000 pixels: 

1 minute of computer time; 

- classification of test zone data: 40 minutes of computer time 
for a real time of about 6 hours of display 4, console in the 
time sharing mode; 

- choice of a printout legend: 2 minutes of computer time; 

- transfer of the classification data to magnetic tape for use 

in the CIRCE plotter: 2 minutes of computer time. 

3.4.2. Use of the HP-2100 system 

- choice of the stipulated area of about 900,000 pixels: 20 

minutes of computer time; 

- cartography at a 1/100,000 scale on the printer-plotter: 

75 minutes of use time. 

3.4.3. Use of the IBM-370 system /15 


- generation of information on the sequence of plots using the 

classification output from the test zone: 1 minute of 

computer time; 

- plotting of a polychrome map at a 1/50,000 scale: 20 minutes 

of use time. 


18 




WWW* 




m ■ii^wiiiinuii.iL'ili mi. nm mi imw» 




CHAPTER 4 - LANDSAT AND NOAA DATA ON CHAD 


4-1* LANDSAT 1 . Data on Chad . 

The N’Djamena area, capital of Chad, having been reserved as the 
test zone, the data used are those of the following three frames: 

- E 1101 08515, November 1, 1972, 

- E 1173 08513 , January 12, 1973, 

- E 1209 08521, February 17, 1973 . 

These three frames have been selected in the series of those recorded 
by NASA like that appearing in the cumulative list issued by the 
EROS DATA CENTER (Figure 4. A). 


The choice of the dates booked was motivated by the idea of com- 
paring the landscape observed between the end of the rainy season 
(November, 1972) , the end of the cool period (January, 1973), and the 
dry and hot mid-period (February, 1973)- Between these three dates, 
the evolution of the countryside is in effect strongly marked by the 
disappearance of the wet areas, the very great lowering of the ponds 
and puddles, the drying-up of the vegetation and the lowering of the 
river waters. 


The N'Djamena zone (12°ia north, 15°00 east) is situated in the 
upper quarter of the two strips 3 and 4 of the frame. 

4-1.1. The photometric parameters 


4- 1.1.1. Image E 1101 08515, November 1, 1972 

- Geographical coordinates of the center: 11°28 N, l4°35 E. 

- Geographical coordinates of the nadir point: 11°24 N, 

l4°39 E. 



- Solar elevation: 51°. 

- Solar azimuth of the radiometer: 132°. 

4. 1.1.2. Image E 1173 08513, January 12, 1973 

- Geographical coordinates of the center: 11°30 N, 14°40 S. 

- Geographical coordinates of the nadir point: 11°26 N, 

. 14 ° 4 6 33 ; 

- Solar elevation: 42°. 

- Solar azimuth of the radiometer: 134°. 

4. 1.1. 3. Image E 1209 08521, February 17, 1973 

- Geographical coordinates of the center: 11°39 N, 14°33 E. 

- Geographical coordinates of the nadir point: 11° 35 N, 

14° 38 E. 

- Solar elevation: 47°. /I 8 

- Solar azimuth of the radiometer: 122°. 

4.1.2. Processing code for the images 

In the following text, each image will be designated by a code 
permitting the assembly of all the processing which pertains to it. 

The image E1101 of November 1, 1972 is called FOR 1; E 1209 of 
February 17, 1972 is FOR 2; E 1173 of January 12, 1973 is FOR 1.2, 
and the processing of dichronic images FOR 1 and FOR 2. 

2. f - The ^ LAND SAT 2. data on Chad 

The data used are those in the following three frames : 

- E 2259 08384, October 8, 1975; 

- E 2349 08374, January 6, 1976; 

- E 2385 08371, February 11, 1976. 


The choice of the reserved dates was motivated by the idea of 
comparing the landscape observed in 1972 - 73 9 during the year of the 
great Sahel drought, with the same landscape observed under normal 
rainfall conditions . 

Analysis of these data will result in the publication of 
dichronic maps placing in evidence the multiseasonal multiannual 
variations such as those resulting from environmental conditions and 
from human actions. 

The photometric parameters are the following. 

4.2.1. Image E 2259 08384, October 9, 1975 
(coded FOR 4) 

- Geographical coordinates of the center: 11°32 N, l4°4l E. 

- Geographical coordinates of the nadir point: 11°32 N, 

14°45 E. 

- Solar elevation: 53° • 

- Solar azimuth of the radiometer: 116°. 

4.2.2. Image E 2349 08374, January 6, 1976 
(coded FOR 5) 

- Geographical coordinates of the center: 11°32 N, l4°35 E. 

- Geographical coordinates of the nadir point: 11°31 N, 

14°38 E. 

- Solar elevation: 40°. 

- Solar azimuth of the radiometer: 133°. 

4.2.3. Image E 2385 08371, February 11, 19_7_6 QS_ 

(coded FOR 6) 

- Geographical coordinates of the center: 11° 33 N, l4°33 E. 

- Geographical coordinates of the nadir point: 11°30 N, 

l4°35 E. 
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- Solar elevation: 42°. 

- Solar azimuth of the radiometer: 122°. 

4.3. NOAA data on Chad 

4.3.1. Choice of dates 

Chad is overflown by the NOAA meteorological satellite every 
12 hours with recording in the morning (Figure 4.B) . With a view to 
thermographic cartography of Ghady ten images have been recorded on 
the descending node between 0720 and O 815 by : choice on those days 
in the months of June and July, 1977 3 on which the axis of the 
trajectory was located between 6°00 east and 22 °Q 0 east at a record- 
ing latitude of 12°1 north. 

The recording was carried out by the APT (Automatic Picture 
Transmission) station of the ASECNA at N^Djamena on analog tapes. 

The analog tapes digitized at C.EwM.S. at Lannion permitted the 
choice of the following dates: 


- orbit NOAA 5 - 3857 3 June 6 , 1977: 

descending node: 7 h-43’26 n , 13 ° 98 E, 

recording: 12°1 N, 17°6 E. 

- orbit NOAA 5 - 3882, June 8 , 1977: 

descending node: 8 h ll , 52 tT , 6° 8 8 E, 

. recording: 12° 1 N, 10°5 E. 

- orbit NOAA 5 - 3894, June 9 . 1977: 

descending note: 7 h 27 r 54", 17°88 E, 

recording: 12°1 N, 21° 5 E- 

- orbit NOAA 5 - 3919, June 11, 1977: 

descending node: 7 h 56 T 193 [sic], 10°78 E, 

recording: 12°1 N, l4°40 E. 
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- orbit NOAA 5 - 3956, June 14, 1977 : 

descending node: 7 h 40 , 45 ,f > l4°67 E, 

recording: 12°01 N 3 l8°30 E. 

- orbit NOAA 5 - 3981, June 16, 1977: 

descending node: 8 h 09 ! 09 " 3 7° 57 Ey 

recording.:. 12" 1 N, 11°2.0 E " 'f-' *' :■>’ 

4.3.2. Characteristic s of the NOAA data /21 

The distance between N'Djamena and the Equator- .along a projec- 
tion of the orbit on the ground being about 1,387 km, the difference 
between the time of passage over the nadir and the time of the 
descending node is about 4 1 (mean velocity at an altitude of 15.10 km 
is 7* 5 km/sec) . 

The analog tape recorded at the N'Djamena APT station (ASECNA) 
was digitized at the C.EJXLS. under the following conditions: 

• Hate of sampling: 1000 points per second. 

•. Number of points per line in the scan: 1250 points; the useful 

portion of the signal corresponding to the effective scan at 
the Earth's surface involves 625 data points among the 1250 
points corresponding to a complete revolution, of the Vraddo-A 
meter mirror. These 625 points cover a scan width of about 
150.0 km. 

• The reolution of the digitized image in the morning is 
7.54 km along the scan and 7.32 lan along the trajectory 
(course of the satellite). 

*' The period of rotation of the satellite is 116 minutes, 

20 second. 


* The pick-off channels correspond to the following spectral 
-bands-: 
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Figure 4. c. Image produced by visible bands (0 
5 satellite. The data received at N'Djamena on 
the north part of Central Africa. Lake Chad api 
m the lower part of the i mag e 
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Figure 4.D. Image produced by the 
pm) by the NOAA 5 satellite. This 
Lake Chad and that of N’Djamena 


thermal infrared band (10.5 - 12.5 
enlargement shows the region of 
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Pi' 


— in the visible: 0.52 to 0.73 microns almost retrieving 

channels MSS 4 and MSS 5 on LANDS AT 3 

- in the infrared: 10.5 to 12.5 microns. 


We show in Figures 4.C and 4.D examples of the printout of the 
NOAA 5 data on the printer-plotter at CTAMN. 
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CHAPTER 5 - AUTOMATIC CLASSIFICATION PROCESSING WITH 
FOR I, FOE 2, FOR 3 AND FOR 1.2 


5.1. Utilization of the FRACAM system 


The classification system developed by CTAMN provides the fol- 
lowing functions: 


- spectral display and statistics (FRAC 1), 

- homologous analysis (FRAC 2), 

- classification in ascending order (FRAC 3) s 

- combining around changing features . 


FRAC 1 provides a general function of display on a TEKTRONIX 
display console allowing the edition of a certain number of images 
and statistical results so as to facilitate the dialogue between 
users of the classification system. This step is particularly im- 
portant at the time of the determination of the cutting out of the 
histograms adapted to the study themes, as well as at the moment of 
the spectral coassociation with the classification system. This last 
phase thus allows one to decide on the pursuit of the cutting out of 
spectral classes, the taxonomic significance of which appears very 
confusing to the user. 


FRAC 2 allows, by means of homologous analysis , the determina- 
tion of a sub-space for representing the binary data associated with 
the cutting out of the histograms . This program follows from an 
ascendant hierarchical classification of slices across the channels 
(FRAC 3) in which the role is to show up thr correlations able to 
exist between the channels at the level of the most significant 
spectral signatures. The Interactive interruption of the classifica- 
tion tree determines the number of spectral classes and decides on 
their characterization by the slices of the channels which, are the 
most c orre 1 at i v e wit h them. 
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FRAC 4 develops the final phase of the classification by 
association of each pixel to the spectral class to which it is the 
closest. The results are then displayed and critiqued by the renewed 
use of FRAC 1, thus permitting 'the] interpretation, of each spectral 
class as a function of the geographic classes on the display console. 

This last step of the geographic and taxonomic post-association 
allows the conduct of the refinement of each of the groups as a func- 
tion of the theme being treated. This approach,, called classifica- 
tion stepping is carried out by reiteration of FRAC 1 - FRAC 4 all 
together in adapting, the statistical parameters used to taxonomic 
.sub-groups which the user wants to discriminate. 

5.2. Data coding 

In order to provide for a degree of acuity sufficient for the 
choice of integrations among the spectral bands, the processing takes 

place by a transformation of the measurements recorded by the MSS 
radiometer. 

This transformation, called disjunctive encoding, is fully con- 
nected with the technique of photoprocessing by equidensities . In 
effect, such a coding is carried out by associating with each physical 
value a binary code corresponding to'the portion of the values (or 
level of gray scale) in which it is included. 

Thus, with each slide of the histogram, there is associated a 
image which, in the methods of colored equidensities, corre- 
sponds to a particular color of the pallette being employed. With 
each slice of the order i, there is associated a relative frequency P^; 

P. = surface covered by the portion of the values i 
1 total surface of the image 

The choice of the cutting is obviously the preponderant control para- 
meter as far as the effectiveness of the corresponding coding is 
concerned. If, a priori, none of the information is available in the 
cartographic area, the choice will perhaps be guided by the information 


30 


theory of Shannon. In this setting, one can characterize the quan- 
tity of information of a disjunctive. coding by the magnitude H: 

•- ' 

: H * - SH Pi 1°S2 p i 
i«l 

where N - number of groups cut. This criterion is maximum when all 
the are equal to 1/N or where, to put it the same way, all the 
binary images are the same surface. This is the case when the cutting 
takes place at groups of isopopulation ; one would then get: 


H 


max 


Iog2 N 


where N = number of groups cut. The quantity H is not, however, a 
mean criterion which is maximized in the case of cutting into iso- 
populations, despite the quite particular significance associated 
with a group of particular value. 


Thus it is that the DUERB program offers, from the selection f 

of data, a cutting into, isopopulation of each of the histograms f 

corresponding to the spectral bands MSS 4 to- MSS 7. Thus, for the 
area FOR 3 (Figure 5. A), the list of the parameters for selection is I 

followed by a graphic presentation of each histogram (Figure 5.B), ] 

accompanied by an appropriate slicing. I 

| ' 

Indivi aually at the level of each binary image of frequency p^, } t 

the preceding criterion can appear to be of little significance. ] 

Association with the spectral quality can thus be practiced in adapt- ! 

the .cutting to the level of taxonomy desired. If one looks at /28 j 

the group of histograms corresponding to FOR 3 (Figures 5.C to 5.F), f ! 

the coding In isopopulation slices can appear acceptable for the ! * 

MSS 4 and MSS 5 spectral bands. On the other hand, the pronounced f 

multimodality of the MSS 7 histogram leads to an adjustment of the 5 

cutting. j 



• Gutting of channel MSS 7 : 
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TAPE IDENTIFICATION 


FRAME IDENTIFIER 

TAPE NUMBER 3 

LENGTH RECORDED 

NUMBER OF POINTS PER LINE 

ADJUSTED LINE LENGTH 


1173 - 0 * 515.00 

*296 

810 

3240 - 


DATE OF FRAME ACQUISITION 12 JAN 73 
COORDINATES OF PHOTO CENTE R 

LATITUDE \C ‘il l -30 
LONGITUDE \ / *01 4 -40 


a. LANDS AT characteristics in Zone 3 


ORDER OF FIRST POINT TO PROCESS BY LINE £ TPDF .9 
ORDER OF LAST POINT TO PROCESS BY LINE . 

NUMBER OF POINTS OMITTED ' N R S A U f - . 0 


280 
5 39 


START LINE " j-OPf 4 ^ 
END LINE l p I v ” 
NUMBER OF LINES OMITTED 


. 1 9 . 3 . 
.34 3 

NISAUT - I' ■ 


NUMBER OF ISOPOPULATION LEVELS N‘»NJV - 


NUMBER OF POINTS READ 'IPTS - 3o2A{) 

" " " " " """" ■"•• tv '- 

NUMBER OF POINTS PER LINE ; I* - 2 AO 


b. FOR 3 parameters if or data selection 

Figure 5. A. List of parameters of the DUERB program for the FOR 3 
area 


The adjustment is used for successive visualization of each 
binary image on the display console, so that the taxonomical signi* 
ficance can be evaluated. 


5.3- Homologous analysis of the image 

The preceding coding allows the conversion of the initial set 
of data into a binary set in which one studies the structure with 
respect to the similarity which can exist between the binary images 
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Figure 5.B. Histogram of channel MSS 4 for FOR 3 area. The values 
in the boxes correspond to levels of cutting out of isopopulations 
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Homologous analysis offers us plane plots in which each binary 
image is represented by a point. Two binary images are all the more 
alike, the closer together the corresponding points are. Each image 
(or each. point) is tagged with a letter/number code made up of one 
letter and one number. The letters from A to D correspond to the 
spectral band going from MSS4 to MSS 7- As for the number, it indi- 
cates the order of the value slice associated with the binary image. j 

The codes corresponding to FOR 3 are indicated in Table 5.G. j 

5 • 4 . Determination of the groups j 


of 


th? 


The phase is that of the homologous analysis and the printout 
binary images . 
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The determination of the groups which are most characteristic 
from the standpoint of their taxonomic significance takes place in 
two stages. The classification tree of images in hierarchical ascen- 
dancy (Figure 5-H) is cut at a constant diameter and one thus visual- 
izes the groups which are formed (the visualization is at the 
1/100,000 scale on the display console, then out to the printer- 
plotter). Then, the study of the parameters for control of the assembly 


TABLE 5-G- SLICES OF VALUES CORRESPONDING TO THE NUMBER-LETTER CODING 
USED IN THE FOR 3 ANALYSIS 














































Figure 5*K. Ascendant hierarchical classification tree for 
binary images 

of and the identification of known objects in the different classi- 
fications allows one to go by successive steps from the uncontrolled 
combinations to controlled combinations by means of spectral support 
(identification of the spectra of known objects) and taxonomic sup- 
port. One ends up with a classification tree and with groups a the 
taxonomic significance of which is certain. 

The process of descending classification has been retained. 

It permits also a by successive iterations a procession to the cutting 
into groups as removed as possible from the function of the recog- 
nition of different objects characterizing the taxonomy. 

From division to division by utilizing the spectra used on the 
statistical parameters of the signatures, one arrives at a legend 
having the most discriminating geographic significance possible within 
the population of pixels (elementary slices of LANDSAT imagery) in 
the zone studied . 
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the control parameters of the assembly: 

example in FOR 3: 

class '91 - ; A1, Bl ' 
class 92 - Cl, Bl, 

a group/ channel cross reference table which one can read out 
on the printer-plotter; 

the spectral statistics for the groups: 


example in FOR 3 : 
class 91 



- graphic readout of the groups in the form of scale-less 

automatic cartography according to the graphic symbols chosen 
in the print channel of the computer. 

Graphic readout and output from the printer-plotter ‘are possible 
at each of these stages. 



CHAPTER 6 . PROBLEMS OP IMAGE GEOMETRY 


6 * 1 * Geometric corrections of LANDSAT images 
(EROS format) 


LANDSAT images are contaminated by a certain number of geometri- /35 
cal distortions, of which the most important is caused by the rota- " ~ 

tion of the Earth during the period of time in which the frame is 
being received (a- 28 sec). 

This phenomenon gives rise to two types of distortion: 

a geometric displacement of lines from the right to the left; 

- a stretching out of the image in the direction from top to 
bottom. 


After analysis, it appears that one can ignore the other types 
of deformation (curvature of the sweep lines, enlargement of the 
Pixels which are at a distance from the nadir point), their magnitude 
m most cases being of an order equal to the pixel itself. 

Line displacement (deformation of the 
parallelogram) 


The fact that the Earth turns during the acquisition of an image 
is troublesome to the extent that the system of synchronization in 
which we have to work is not fixed. Figure 6. A shows by way of 
example the distance BC, giving the maximum distortion due to the 
rotation of the Earth during the 28 sec of the image acquisition or 
0.11 degrees, which gives 13 km at the Equator. 


A Canadian publication* allows us to formulate the problem. 


Cartographic accuracy of ERTS. Dr. 
Council of Canada , Ottawa, Canada. 


V. Kratky. 
KIAOSI. 


National Research 
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One finds 


w ■ C \l sin2 i - sin^^.. 

where i = orbital inclination of the satellite (i=99); 
cj) = latitude of point A; 

C = constant (1/13-9) connected with the speed of movement of 
the satellite. 



Figure 6. A. Representation of the deformation of a parallelogram due 
to the rotation of the Earth 

6.1.2. Stretching out the image in latitude 

This distortion is the direct result of the deformation of the 
parallelogram . 


If the Earth had not moved during image acquisition, point B 
would get position C (Figure 6. A). Under these conditions, the image 
stretch out is the amplitude DC. 

One finds: Sc =* AC x G x sin (i -tf/s) • 

AC is the coarse height of the image (185 km), and C is the 
constant defined above (1/13-9). It is noted that this quantity DC 
is a constant value = 2.04 Ion. 

6.1.3. Generation of the image on the VERSATEC printer 

On the VERSATEC printer, each pixel is represented by a rectan- 
gle of 4 x 6 VERSATEC units (i.e., 0.5 x 0.75 mm) for a scale of 
1/100,000, and 8 x 12 units (i.e., 1 x 1.5 mm) for the scale of 
1/50,000. 

Under these conditions, it is necessary to begin with to have 
an image reformatted such as to get pixels representative of a sur- 
face of 50 x 75 meters. This work is performed at the same time as 
the distortion correction analysis above. 

6.2. An experiment in marking the LANDS AT image 
with a mirror 8 


An experimental mirror (Plate B) has been installed at the 
Science Institute in N’Djamena 88 (April, 1977). Two sightings were 


The experiment is based on that of W. E. Evans, described in 
"Marking LANDSAT Images with Small Mirrors Reflectors" in NASA 
Earth Resources Survey Symposium, Houston, Texas, June 1975, Vol. 
I.B, pp. 1185-1196. 

## 

The experiment was carried out in collaboration with J. P. David, 
Professor of Chemistry, Dean of the Science Faculty the Univer- 
sity of Chad, and Leandri, maintenance engineer at the Science 
Institute. The position of the mirror was directed by G. Doublier, 
geometrist and Professor of the E.N.T.P. at N’Djamena. 
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simulated on May 24, 1977 and June 11, 1977 , the installation being 
provided with a temporary mirror. 

6.2.1. The observation of May 24, 1977 

The data from the simulation are the following: 

- assumed time of LANDSAT satellite overpass: 8 h 37 f 04” UT; 

- assumed altitude of the satellite: 912.4 km; 

- mirror-satellite azimuth at the moment of the sweep: 304 gr 

26692 ; 

- azimuth of the sun at 8 h 37' 04": 78 gr 7453* 

6.2.2. The observation 
of June 11, 1977 




The data from the simulation 
are the following: 

- assumed time of LANDSAT 
satellite overpass: 

8 h 37 T 04” UT; 

- assumed altitude of the 

satellite: 912.4 km; 

- azimuth of the bisector 
of the satellite-mirror- 
Sun angle with respect to 
geographic north: 

394 gr 4126 

- azimuth of the Sun at 
8 h 37' 04" UT: 

74 gr 5582 

The mirror was installed 
according to the following co- 
ordinates : 


Mirror for marking LANDSAT images — 
May, 1977 

Plate B. Faculty of Science, 
University of Chad 

(Plate continued on following 
page) 
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Positioning of the mirror 
Plate B. (continued) 


• latitude N: 12° 07 1 41', 

• longitude E: 15° 00’ 11”. 




I 




6.2.3. Problems posed by the actual observations 

To arrive at useful results, the following problems need to be 
resolved : 

a) knowledge of the ephemerides of the LANDSAT satellite; 

b) knowledge of the precise altitude of the satellite (varies 
between 909 and 920 km); 

c) acquisition of an aluminum mirror especially polished for 
the best possible reflectance. The mirror must be slightly 
convex ; 

d) the availability of atmospheric refraction measurements and 
of atmospheric pressure. 
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CHAPTER 7- GEOGRAPHICAL TAXONOMY 


7-1* Ground truth documentation 


/40 


Geographical features were made the object of an inventory in an 
area 100 km around N’Djamena for which the maximum number of docu- 
ments were brought together. In particular, they comprised: 

7-1-1. Basic topographic maps of the I.G.N.F, 

The area studied is covered by the following maps: 

- 1/1,000,000 sheets of Fort-Lamy NC 33 and of Garoua ND 33, 

- 1/500. COO sheets Lake Chad (Maine-Soroa) and Fort-Lamy, 

- 1/200,000 sheets Fort-Lamy ND 33 IV, Makari-Goulfey ND 33 III, 
Mora NC 33 XXI, Mogroum NC 33 XXII. 

7-1-2. Aerial photography of the I.G.N.F. 

The area studied is covered by the following flights: 

- AE 391/500 - 1971 - March 23, 7 h to 12 h; 

- TCH 02/150 - 1972 - March 19, 10 h to 11 h; 

- TCH 74 08/200 - 1974 - July 11, 12, 13; 9 h - 12 h; 

- TCH 74-75 10/500 - 1974 - 1975 - Mogroum sheet, photos 
5929 to 6665 on 29/12/74 to 10/1/75, 9 h to 12 h. 

7-1-3- Land registry plats 

They are of a scale of 1/2,000, 1/5,000, 1/10,000, and 1/20,000. 
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The area studied, is covered by the 1/200,000 sheets of J. Pias, 
E. Guiehard, and J . Barbery, dating from 1963. They are the sheets 
of For t -Lamy, Port -Four eau, Mora, ; and Mogroum. 


7 . 1 . 5 . 


geological map of the D.M.G of the A.E.F. 




i 

This is an exploration map at 1/1,000,000 laid out by J. Barbeau | 

in 1956 for the D.M.G. of the A.E.F. ; 


7.1.6. Soil data 


C. Tobias focused on these data in March, 1976, in technical 
note no. 13, O.R.S.T.O.M. - N'Djamena, "Soil Utilization and Irriga- 
tion Capabilities in the N'Djamena Area", 14 pages, mimeo. 


7 . 1 . 7 . 


Surveys based on ground information and 




These have been put together by C. Bardinet at the Laboratory 

for Cartography of the University of Chad from 1973 to 1977. j 

{ 

j 

t 

t 

Three types of documents have been used on the ground: i 


• aerial photographs in I.F..C. and in panchro (N'Djamena 
Aeroclub) and ground photos in I.R., I.R.C., and panchro 
taken by C. Bardinet from 1973 to 1977; 

• the photointerpretation surveys after the aerial coverage of 
the I.G.N.F.; 

• photointerpretation, surveys following the enlargements of 
LANDSAT image E 2349 08374 of January 6 S 1976, enlargements 
carried out by Gaudin and Husberg at the C.E.R.C.G. of the 
CNRS using the 1/1,000,000 negative from NASA. These enlarge- 


ments were prepared, channel by channel, at 1/200,000 and at 
1/500,000 for the N'Djamena area. A composite color 




photo has likewise been prepared on transparent film at 
1/200,000 on Kodalith film (Plates I, J, K, and L). 

7.2. Locating; the taxonomic subjects 

A certain number of features have a pure signature, for example, 
river waters, bare and dry sands. Other features have a mixed sig- 
nature, for example, vegetation which may be more or less chloro- 
phyllic, more or less associated with water or with very humid areas, 
or again the urban areas in the case of which the zonal division into 
quarters can be marked either by the presence of green areas in the 
modern residential quarters with villas and gardens, or an earthen 
habitat (the poto-poto) more or less crowded and with little or no 
vegetation. 

In the test area of around 200 x 200 pixels, the features have 
been located, and with each feature a certain number of objects have 
been specifically indentified, the spectral signature of which, by 
their development through dichronic analysis, should let us refine 
the support of the classification process. 

7.2.1. Pixels, features, and objects 

On the LANDSAT digital tape, the features and the objects which 
they compose are recognized by the form of the levels of mean lumi- 
nance which is recorded in the 79 m x 79 m pixels, which with overlap 
would be 79 m x 59 m. The index of the pixels is selected in channel 
MSS 5 from the photo which is delivered with each tape. 

7 . 2 . 1 . 1 . The index of pixels FOR 1, FOR 2, FOR 3 

• The FOR 1 index (November 1, 1972, strip 3) 

lines 100 to 250 , 

1 st pixel: 401, last pixel: 65 9 > 

number of pixels read: 39 , 109 , 

number of pixels on a line: 259 . 
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photo has likewise . been prepared on transparent film at 
1/200,000 on Kodalith film (Plates I, J, K , and L) . 

7*2. Locating the taxonomic subjects 

A certain number of features have a pure signature, for example, 
river waters, bare and dry sands. Other features have a mixed sig- 
nature, for example, vegetation which may be more or less chloro- 
phyllie, more or less associated with water or with very humid areas, 
or again the urban areas in the case of which the zonal division into 
quarters can be marked either by the presence of green areas in the 
modern residential quarters with villas and gardens, or an earthen 
habitat (the poto-poto) more or less crowded and with little or no 
vegetation. 

In the test area of around 200 x 200 pixels, the features have 
been located, and with each feature a certain number of objects have 
been specifically indentified, the spectral signature of which, by 
their development through dichronic analysis, should let us refine 
the support of the classification process. 


7-2.1. Pixels, features, and objects 

On the LANDSAT digital tape, the features and the objects which 
they compose are recognized by the form of the levels of mean lumi- 
nance which is recorded in the 79 m x 79 m pixels, which with overlap 
would be 79 m x 59 m. The index of the pixels is selected in channel 
MSS 5 from the photo which is delivered with each tape. 

7-2. 1.1. The index of pixels TOR I , FOR 2, FOR 3 
The FOR 1 index (November 1, 1972, strip 3) 
lines 100 to 250 , 

1 st pixel: 401, last pixel: 659, 

number of pixels read; 39,109, 
number of pixels on a line: 259 . 
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PLATE J N'DJAMENA AREA 1 

Enlargement of the EROS CENTER 1/1,000,000 negative by MM. Gaudin and Husberg at the photo 
lab of the C.E.R.C.G. of C.N.R.S. in October 1977. LAND SAT image 2348 08374 of 1/6/1976 

Enlargement to 1/200.000e Photo 2 - channel MSS7 - 800 - 1100 micron 









































negative by MM, Gaudin and Husherg at the photo 
ir 1977. iiNDSAT image 2348 ;d83;74 of ;l/6/1976 








































■ The FOR 2 index , (February 17 1973, str .p 3 ) 


lines 354 to 


number of pi 
number of pi 


i last pi 
read: 39 * 1 ^, 

on a line: 259 . 


* The F.QR 3 index (January 12, 1973, strip 3) 

lines 193 to 343, , / 

1 st pixel : 28,0 jy last pixel :■■■ 538' , 

number of pixels read: 39 , 109 , 

number of pixels on a line-: 259 . 


7 . JPhe dichronic index for pixels in FOR 1.2 . 

An index of 150 lines and 259 columns has been created in which 
one takes into consideration the diehr onic signature of the pixels of 
FOR- 1 and of FOR 2. The line n and -tie column' m for FOR 1 are super- 
imposed on the line n+25,4 and column m +151 of FOR 2.. This super- 
imposition is possible when all of the indices have received the cor- 
rections necessary for them to present the same geometry. 

. vv'f: .• , 712.2. i Feature a, andv -refer one e vob 1 ects 


ft 


7 • 2 . 2 . , 1 . Alignment of pixel ; reference 


The width of each basic line is 79 m on the ground and the 
information collected in each spectral band corresponds on the ground 
to . a square of 79 m on a side. The dimension of this basic measuring 
square is prescribed with account taken of the orbital altitude of the 
satellite by the instantaneous optical field which is 0.086 milli- 
radians. .This optical field which is square in form has. the same 
dimensions for each of the spectral bands MSS4 to MSS 7- 
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The frequency of oscillation of the spin-scan mirror (13.62 Hz) 
combined with the nominal satellite orbital' velocity (6.47 km/sec) 
assures an almost perfect contiguity of the basic lines at a 79 m 
width. 

This contiguity is not maintained in the direction of the scan. 
In effect , the signal, having been sampled at 9-95 microseconds, the 
distance between the measurement sampling points is 56 m. 


Prom this one can infer that the pixels or resolution elements 
sketch out a regular rectilinear grid (Figure 7. A) on the ground 79 
m x 56 m, in which each node is associated with an elementary radi- 
ance measurement corresponding to a square of 79 m on a side. 


79 m 



+ 4 - + + 


Figure 7. A. Regular rectilinear pixel grid and basic 
measurement square. 

A surface of 1 km 2 is covered by about 224.28 pixels because /48 

there are 12.6 pixels per kilometer in the direction of the columns 
(velocity vector of the satellite) and 17.8 pixels per kilometer in 
the direction of the scan lines. 
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7 . 2 . 2 . 2 . Effective spatial resolution of LANDSAT 


To choose the objects marked for identification and locating 
on the ground, it is important to take into account the limit of the 
resolution of the MSS radiometer. First of all, it is appropriate 
to separate two distinct problems which the experimenter can en- 
counter simultaneously: 

• the capacity of the MSS radiometer to distinguish two objects 
of the same type which are located on the ground at a certain 
distance from each other is essentially concerned with a 
problem of detection ; 

• the way in which the MSS radiometer follows the spectral con- 
trasts of the terrain observed and reproduces with fidelity 
the spectral signature of the objects for the cartographer. 

The resolution limits corresponding to each of these problems 
are quite different. Therefore, we propose to clarify here the con- 
cepts which apply to them. 

a - Limit of resolution in detection 

Conforming to the ERTS Data User Handbook No. 71SD4249, we will 
confine ourselves to the problem of the detection of a binary terrain 
feature consisting of three parallel rectangular strips (Figure 7.B-1) 
following the model of the USAF tri-bar target, 1951. In order that 
we put ourselves in as realistic conditions as possible, we have 
chosen terrain features (Figure 7.B-2) consisting in part of objects 
quite often encountered in the LANCHAD project and enabling an evalu- 
ation of the possibilities of detection with the MSS in situations 
of vegetation/bare earth contrast. 

From the point of view of the vegetation, we have chosen both a 
shrubby carpet of vegetation common in the savanna area and a seeded 
area of grain field having a growth presenting the mean proportion 
of vegetative covering corresponding to the cultivation means 
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and N’^ are the 

reflectivities of two 
types of objects 


Figure 7.B-1. Model of the binary terrain feature 
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7.B-2. The four types of terrain features considered 

Figure 7.B. Theoretical model used for the estima- 
tion of the detection resolution limit of the MSS 
radiometer 
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generally employed in Chad. For the soils, we have retained a sandy 
soil of strong reflectivity and a wet argilaceous soil of dark color 
able to represent the areas of alluvial cultivation. The reflec- 
tivity spectra of these four objects are presented in Figure 7.C. 


In the model presented, a certain number of parameters, depend- /52 
ing on the MSS spectral band used, should be considered in determin- 
ing the limit of the detection resolution; the question is: 


- (MSS^) = the reflectivity of the most reflective object in 

the MSS i spectral band, 

- N 2 (MS3 i ) = the reflectivity of the least reflective object in 

the MSSi spectral band. 


- C. 
1 


" N. 
1 


N 2 (MSSi) 

V' (MSS ) 

1 1 band , 

N 2 (M33 i ) + N x (MSS.) 


= reflectivity contrast for the MSSj_ spectral 


mean reflectivity 


of the terrain feature in the MSSj_ spectral band 


Tables 7.D and 7.E present the values and the different para- 
meters for the terrain features obtained, which already enables one 
to make a certain number of statements concerning the variation of 
contrasts as a function of the spectral band. One notes that no 
distinction is possible between crops and argilaceous soils in MSS 5 
(0.6 - 0.7 ym) , but that, on the other hand, the best contrast is 


TABLE 7.D. TABLE OF THEORETICAL REFLECT- 
ANCES CORRESPONDING TO THE SPECTRAL BANDS 
OF THE MSS RADIOMETER 



MSS 4 

MSS 5 

MSS 6 

MSS 7 

Sand 

25 

34 

41 

46 

Black argile 1 

5 

7 

10 

15 

Cultivation 

9 

7 

23 

25 

Shrubbery 

9 

12 

25 

35 
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TABLE 7.E. TABLE OP CONTRASTS (C) AND OP AVERAGE REFLECTIVITIES (N) 
CORRESPONDING TO A REGULAR ALTERNATION OP TWO TERRAIN FEATURES. 

THE CONTRASTS MARKED BY A RECTANGLE RELATE TO THE CASE IN WHICH 
THE SOIL IS MORE REFLECTIVE THAN THE CROPS 



MSS 4 

MSS 5 

MSS 6 

MSS 7 


C 

¥ 

m 

N 

C 

N 

C 

S 

Shrubbery/ 
sandy soil ; 

2,78 

17 

2,83 

23- 

1,64 

33 

1,31 

40,5 

AS 

Shrubbery/ 
clay soil ; 

1,80 

H 

1.71 

9,5 

2,5 

17,5 

2,33 

25 

AA 

Cultivation/ 
sandy soil 

2,78 

17 

4,86 

20,5 

1,78 

32 

1,84 

35,5 

CS 

Cultivation 
clay soil 

1,80 

B 

1,00 

a 

2,3 

16,5 

1,67 

20 

CA 


obtained in this same band between shrubbery or crops and sandy soil. 
Concerning the contrast between shrubbery or crops and argilaceous 
soils it is the MSS 7 channel (0.8 to 1.1 pm) that one should choose. 
Finally, the exact contrast of a type of an object with respect to a 
variable natural environment depends greatly on the spectral band 
employed . 

In accordance with the ERTS Data Users Handbook, we will define 
the detection resolution limit in the width of the strip &, expressed 
in meters (Figure 7-B-l). The latter at the same time depends on the 
contrast (C) and the mean reflectivity (N) , and therefore also on 
the spectral band used. This dependence becomes clear by means of 
the transfer function (modulation transfer function) and the signal 
to noise ratio, two mathematical concepts which we will not explain 
here. However, it is possible to work up graphs giving Z as a func- 
tion of the different parameters C, N, and MSS ± . We present here 
two examples of this type of graph. 

In the first example (Figure 7.F), the contrast has been fixed 
at C = 1.8; it is the case of features CS in MSS 7 and of CA or AA 
in MSS 4. One gets: 
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-& CA - = 47 m in MSS 4, 

4 cg 49 m in MSS 5- 

These lower resolution limits at 79 m may surprise the reader. 

It is connected in ;f apti with the optimum capabilities of the system 
concerning which information theory tells us that the lower limit 
of resolution is located at % * 4-9-5 m, even though the width of 
the measuring unit is 79 m. 

The sampling rate of the signal is a supplementary factor which 
practically limits the detection resolution to 56 m in the direction 
of the lines and 79 m in the direction of the columns on the raw data 
tapes received from the EROS DATA CENTER. However, a theoretical 
lower limit at one of the two figures given above is an incentive 
to re-sample the data by interpolating between the pixels in order 
to facilitate the detection. 

In the second example (Figure 7.G), the aim is to show, by means 
of a second spectral band (MSS 5 5 for example), the combined effect 
of the contrast and the mean reflectivity at the resolution limit. 

We already know that in this case the feature CA appears as a 
uniform field (&q A “ 00 in MSS 5)5 in the other cases: 

- *= 45 m in MSS 5 a 

- i Ag = 46,5 m in MSS 5 , 

- 1, = 61 i in MSS 5. 

AA 

On the other hand, when & is larger than these, re-sampling 
will not result in any supplementary information. 

b — Limit of resolution in spectral fidelity 

In most of the methods of data processing, in remote sensing, 
the experimenter must conduct a spectral support phase in the course 
of which a double unequivocal relationship must be created between 












the type of object and the spectral signature which characterizes it. 
To the problem of the detection of objects, there is therefore 

joined that of the ^purity” of the spectral signatures thus 
collected. 

We a sain engage ourselves in the matter of the binary features 
described above in asking ourselves which circumstances will assure 
that the spectral signature of a cultivated plot from the MSS radio- 
meter is most characteristic of this type of object. 


The response-.to this problem is given by a curve called square 
Kave_response 3 giving the maximum amplitude of the signal as~a func- 
tion of the value L of the rectangular bands of the feature model 
adopted. This curve, independent of the wavelength of the observa- 
tion, is presented in Figure 7 -H. The absolute maximum amplitude N 1 
observed in the output of the MSS radiometer is obtained by a value L 
determined by multiplication of the actual mean reflectivity N by a 
corrective term a(l) £ 1 ; 

N* 1 - A (L) x N 

There will then exist a diminution of the spectral fidelity of 
the MSS radiometer as a function of the shape of the object observed. 
This phenomena is illustrated in Figure 7. I where, for an input sig- 
nal corresponding to a succession of two types of objects marked by 
the respective reflectivity N 1 ^ and N ! 23 the maximum (amplitude of the 
radiometer measuring signal varies as a function of the shape of 
these objects. 

According to the characteristics and the performance of the 
classification system used, the spectral support phase imposes a 
fidelity threshold below which the spectral signature cannot be con- 
sidered pure. In general, one will select a threshold A = 80 $, 
which sets the spectral fidelity resolution limit L at 100 m . 

In conclusion, it is necessary to distinguish, at the level of 
the utilization of LANDS AT data, the resolution limits for detection 
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Max. relative amplitude, % 



Figure 7-H. Square wave response of the MSS radiometer 


and those for spectral fidelity. The first come into play in the 
definition of the geometry of objects and can give rise to the 
interpolation of the data in certain cases. The second impose a 
limit on this side of which one can come down without strongly 
changing the quality of the spectral information. The most unfavor- 
able case is encountered when the detection resolution is worse than 
that of the spectral fidelity (A < 100 m) , as for example is the 
case in MSS 7 with thickets of shrubby vegetation on argilaceous- 
sandy soil. Finally it should be noted that the values presented 
in our paper correspond to a theoretical mode; they must be con- 
sidered to be the upper limit of the performance that one can expect 
from LANDSAT. 
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Figure 7.1. Diminution of the spectral fidelity of 
MSS as a function of the size of the objects as 
observed 

The solid line is the input signal; the dotted line 
is the MSS response 



The principal features of the test area are the following 


- river waters, 

- standing water in ponds and puddles, 

- flooded areas, 

- green vegetation, 

- bare sand, 

- urban quarters of traditional poto-poto habitations, 

- urban quarters of modern habitations, 

- cultivated areas on burnt-over terrain. 
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This general taxonomy of the area enables, identification of the 
test area and guides the simplest choice of the classification. 


b - The recognition and ‘the location of features 


The recognition and location: of features must be the more pre- 
cise seeing that one wishes to carry out the dichotomy of a group, 
for example, that of very, green vegetation in which wooded areas 
and forests, are distinguished, on the one hand,, and meadow areas or 
bottom land cult i^^ion^'.©^ 


In the successive steps of classification, each feature must be 
recognized with the maximum precision allowing the interpretations 
of the variations in spectral statistics. 


This very precise taxonomy allows two types of support: 


an internal support beginning with the recognition of the 
spectral characteristics of the processed data; this permits 
decisions about the choice of classification; 


an external support enabling one to take account of other data 
than that which one is processing, notably the knowledge ac- 
quired in the investigations of ground truth concerning the 
possible or probable evolution of the features during the 
dichronic analysis, i.e„. 


• change from water or wet areas into dry, 

• the agricultural calendar, etc. 


Geographic support allows both the precise determination of the 
limits of a feature studied and the forms of the seasonal evolution 


of this feature. 


c - N^Djamena area 




In this area one has determined the following geographic features: 
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- sands: • bare and dry, 

* * — 

• dry with vegetation .of low chlorophyll 
e '•damp , 

• wet ; 

- waters: • flowing rivers of many degrees of turbidity 5 

• standing water of ponds., puddles, and flooded areas; 

- water and. vegetation complexes ('Elates C, D, F) : 

• ponds .with bushy vegetat ion o r with Wat er lilies, 

. • shrubby flooded areas; 

- damp or' wet soils with transient vegetation; 

- wooded areas of high chlorophyll (Plate E); 

- soils dominated by spiney vegetation: 

• black tropical argil es , 

• argiles.with chalky concretions, 

• argile sands, I 

• argile silts, ... . j 

• sandy argiles; j 

- sandy ridges of the bank or of the old hydrographic networks; i 

- cultivated areas, (Plates D and E) : ] 

• dry land cultivations on burned-over areas, j 

• flood plain cultivation, j 

i cultivation in rainy areas i 

I 


• bottomland cultivation 

4 fruit tree cultivation; 

- urban areas (Plates G and H) : 

• areas of modern construction with strong reflectivity, 

• areas of Modern habitations — villas and gardens, 

• landscaped "urban areas of trees with high chlorophyll, 

• bare urban spaces of high reflectivity (stadium, 
racetrack, square . . . ) ■ 
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Plate C. Koundoul Road, October 29, 1975 

1 — panchromatic; 2 — infrared; 
foreground — pond with water lillies; 
middle distance — trees with high chlorophyll; 
rear — copse of acacias 
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Plate D. Koundoul Road, October 29, 1975 

1 — panchromatic; 2 — infrared; foreground — pond with spiney plants and 
water lillies; middle distance — field of manioc; rear — field of sorghum 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 
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Plate E. Koundoul woods (edge of natural second growth) 
October 29 , 1975 


1 j panchromatic; 2 — infrared; foreground — blacktop 
road; middle distance — right, protected replanting: 
left, sorghum field 
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Plate F. Replanting of eucalyptus, Koundoul Road, 
October 29, 1975 

1 — panchromatic; 2 - infrared; foreground — flooded 
area, stalks of sorghum projecting from the water; 
middle distance — replanting of eucalyptus in 1.2 m 
of water; rear — flooded area 


70 



2 


Plate G. N’Djamena, December 14, 1975, site of the 
industrial secondary school 

1 — panchromatic; 2 — infrared; foreground — bare 
soil and dirt trail (argilaceous sand; background 
— the school, aluminum sheet roof, traditional 
poto-poto habitation 
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Plate H. N’Djamena, December 15, 1975 

1 — panchromatic; 2 — infrared; foreground — 
bare argilaceous sand and concrete posts; middle 
distance — observatory dome and hedge of high 
chlorophyll 


7 . 2 . 2 . 4 . Ob j eats identified': in the taxonomy 

i^Djamena area 

In carrying out a descending classification, it Is indispens- 
able to locate the significant objects for each type of feature. 

In the N'Djamena area* the following objects have been located: /66 

- pld meanders and beds of the Chari going from damp to very 
wet with green vegetation (brush or cultivation) , 

- ridges of the bank present and past (sand berms of great 
reflectivity ) 5 

- present and past drainage networks, 

- green areas., woods 3 and gardens at Milesi and Farcha, 

- green urban zones of quarters and boulevards planted with 
trees, 

- urban quarters according to the type of habitation and their 
network, 

- urban spaces of specific reflectivity, 

- agricultural type areas. 

A more precise nomenclature for these objects classified by 
feature can be described as follows: 

Urban features (Plates G, H, and Figure 7-J) 

• urban quarters dominated by greenery: Djamb and N'Gato 3 

Rond-Point, Etoile, Matembe where a habitat of villas with 
gardens and tree-^pdh^avenues predominates; 

• the Bololo, Djamb el Barr quarters of very closely packed 
network, with trees; 

• the Avenue Mobutu 3 University 3 Eboue Secondary School 
district dominated by tree planting; 

- the Gardole, Sabangali quarters where the Italedyl devel- 
opments predominate (strongly reflective roofs and gardens 
with shrubbery); 



Figure 7.J. Map of N’Djamena at 1/100,000 with district toponymy, 1976 


Legend: quarter boundaries; . — quarter names; 

quarter boundaries; ••• — limit of conurbation; +++ 

frontier; principal roads; - — — limits of 

1 " rivers; 9 — central market - mosque 


— sub- 

— international 
the banks; 


A - Gardole 

B — Ambassatna, Kabalal • 

C - Blablim, Centre 

D — MardjSm Dafaq 

E - Bololo Djamb al Bahr 

P - Sabangali 

G — Ardeb DjomSl 

H - Ridina 

1 — Soldiers home ; 

J - T ’iile Arab, Sarro, Hille Leolerc 
K — i. uxa Al Hamra 
L - Kileb MSt 

M — Djanb an N'Gato, Matambe, Etoile 


N — Airport 

0 - Moursal, Paris Congo 
P - Chagoua 
Q - Dembe 

R - Zongo, Grediya (radio stn.) 
S - Am Riguebe, Hille Gazfiz 
T - Cholera Clinie, 40th St. 

U — 30th St. (Gozator) 

V - Diguel 
W — Naga 

X — Syrte 

Y - Farcha, Milgzi 
Z - Ocam 



• the Cite O.C.A.M. district where there is a mixture of 
development, bottomlands; reforestation; 

• gardeners T ponds; 

• the I.Z.V.A.C. gardens and the Farcha cemetery; 

• the airport area; 

• the urban spaces of strong reflectivity , racetrack, St. 
Martin plain, the Gredyia antennas; 

• the central, market district with the great mosque and the 
House of Congress; 

" the Chari bridge at Chagoua; 

• the traditional .residential districts around courtyards, 

arrange in a close network, in particular: Ambassatna, 

Kabalai, central Blablim, Mardjan Dafaq, Ardeb Djomal, 
Ridina, Bout al Hamra, Hille, and Sarro; 

• the traditional residential districts with scattered dwell- 
ings, in particular: Naga, Adiquel, Am Riguebe, Gredyia, 

on the one hand., and, on the other hand, Sara Moursal, 
aris-Congo, and Chagoua west; 

• the nearby uplan districts Adiguel North, Chagoua East, 
Dembe, Goro Waga on the east and, on the west, Farcha and 
Milesi; 

• Kousseri (formerly Fort Foureau) in Cameroon, partly along 
the Chari in the category of modern urban districts with 
tree-lined avenues and partly with traditional residences 
tightly packed in the middle and scattered out to the south. 

Vegetation features with high, chlorophyll : 

• -Milesi woods; 

• the Kouseri woods, and the Kalamaloue Reserve in Cameroon; 

• the Koundoul woods; . . 

• the bottomlands and orchards of Goro Naga, Gassi, and 
Mand Jaffa; 

■ the Digangali bottomlands; 

• modern irrigated farms at Koundoul, Malo in Chfad, and at 
N’Douf in Cameroon. 
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Sandy- features and ridses of banks 


... the Irounda flood area- downstream from Milesi; 

• the Chari flooded areas i in particular near the bridge at 
the confluence with theLagone. • 

Features of the old., meanders of the Chari, at Dignagali and 


Lake and, pond features , in particular at Dembe and NJair, to 
the north and east of N T DJamena, in the close-in urban zone and 
to the southwest of Kousseri and of Mont Ct. Lamy . 

The objects which were located can be classed into three types: 

- those where; the reflectivity varies little over time, 

- those where the reflectivity is complex and variable, 

- those where the reflectivity did not vary during the 
observation period. 

SfeWlSPiBi at 1 / 100 . 000 : geographical taxonom 


of the N'D.jamena- Kousseri area 

This map was created starting with a selection of data in the 
FOR 3 classification described -in section 8 . 3 . 3 . The basic data 
are those of LANDSAT 1, image E 1173 08513 of January 12, 1973- 

The refinement of the controlled classifications allows finally 
the classifying of the basic numerical data as a function of their 
taxonomical significance. This classification can only be created 
on an index of reduced dimension — in this case, 40,000 pixels. 

The question which arises here is that of the extension of the 
classification to all four strips-, of the image under study. This 
question cannot be resolved directly 5 and one cannot consider 
a Pply3-ftg the classification created in the test area to the entire 




7.3.1. Identification- through the classification of 
significant objects in the taxonomy 

Beginning with the. controlled classification of the test zone, 
it is, on the other hand, possible to identify the spectral signa- 
ture of certain key objects in the taxonomy. 

One of the objectives is to jme^£g^|; from the perspective of 
land management : 

- the evolution -of -the- cultivated areas, 

- the extent of the burned-over areas connected with tradi- 
tional agriculture, 

- the spatial extent of beneficial vegetation in the rural and 
urban economies , 

- diminution of the vegetation on account of the intensive 
exploitation, carried out by stockmen, farmers, foresters, and 
urban developers. 

In the four channels, the selection of the reflectivity level 
for these objects is possible starting' with the spectral statistics 
of the various classification groups. This technique makes use of 
electronic filtering means which isolate, in the four channels of 
the image, the objects according. to their specific signature in two 
or more channels. The limits of* the coding can be considered to be 
an upper and a lower filter. 

On the map which is annexed, there have thus been determined: 

- flowing waters exclusive of all lakes and flooded areas which 
can be cultivated during subsidence, 

- burned-over areas which are signs of man's actions in the /70 

countryside — actions by herdsmen and farmers, 

- vegetation areas which one observes according to the more or 
less great intensity of the chlorophyll effect. These allow 

'K^vaiB'the location of dense forestation, diminished forestation and, 
depending on the progress of the vegetative cycle, cultivated 
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zones and pasturage, 

- sandy areas without vegetation which mark the uncultivated 
and arid spaces in the countryside. 

One of the characteristics of this cartography lies in the 
selection of the patches or pixels as a function of a thematic 
analysis. All the pixels not concerned with*- this analysis are, 
in the map presented in the annex, represented by a white space. 

This is to say that, at the maximum., it is possible then to 
create a thematic map considered, useful from the point of view of 
land management so long as the spatial extent of the phenomena 
observed and measured is characterized by the level of specific 
reflectivity. 

• 3 2'-7 Cartography using the satellite image data bank 

Finally, one here approaches the possibility of creating a 
data bank beginning with the indices constituting the recordings 
made with artificial Earth satellites. 

This data bank is only possible if the analysis -through con- 
trolled classification, in the case of a well-located geographical 
region, has allowed us to establish a good relationship between 
spectral statistics and geographic taxonomy. 

But, when this objective has been reached, one can consider 
that for the area studied a geographical data bank is on the way to 
being constituted. This ought to allow one to consider the later 
use of data outside that of the digital image which results from 
the collection of qualitative and quantitative information localized 
to the ground. 

The data bank used for land management appears now to allow a 
thematic cartography resulting from the combined analysis of digital 
data from the satellites and of data collected on the ground by 


direct investigation or by other methods as for example photography 
in black and white* color, infrared,* and false color. 
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CHAPTER 8 - CHRONOLOGY OP COASSISTANCE 
PROCESSING 


■ • 


The descending classification method was used for establishing 
legends of FOR 1, FOR 2, and FOR 3, and then for FOR 1.2, using spec- 
tral coassistanee and post assistance taxonomy. We will summarize., 
the processing sequence for each map by using several selected exam- 
ples. This allows us to explain the problems encountered in the iden- 
tification of the main feature objects as a function of the spectral 
statistics on the one hand, and ,qf< the binary image cartography at the 
scale of 1 :100,00($'and 1:80,000°. 


8.1 FOR 1 Processing Sequence 

The classification in steps, using successive cutting off of 
groups (spectral coassistance) as a function of their taxonomic im- 
portance, is based on the analysis of transformed canal boundaries and 
their interpretation on the binary image. 

The following two tables summarize the classifcation steps using 
five groups, and then other sub-groups. 


8.1.1 Classifications of FOR 1 


Groups Sub Groups 


Transformed Canals 
and Reflection 
Levels 


Features 


I 

A8 B8 C9 D8 Sandy 

66 73 71 31 . 

21 

A1 B1 Cl D1 . Water and Puddle 

40 35 39 18 Vegetation 

22 

C2 D2 A5 B4 Puddle Water 

43 18 47 45 

23 

C4 P4 i Dense Vegetation 

49 24 3 Forests 

. 24 ’ 

B2 C3 D3 Bush Vegetation 

38 45 21 


Cl D1 River 

n Water 

35 ll 

41 

04 D4 A4 B4 Sandy Argiles 

66 30 61 67 
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43 j 

Cl D1 A2 B2 
60 27 58 62 

• Traditional central 
urban area 

44 . 

> • 

A1 B1 C2 B2 
50 60 61 28 ’ 

, Dry Argiles , dry 
vegetation 

51 

A1 Bi' Cl D1 
45 44 50 24 

Green village, modern 
urban area 

52 

C1.B1 B2 
49 22 48 

Black argiles 3 cultures, 
vegetation 

53 ‘ 

A2 B3 C2 D1 ■ 
49 49 53 25 

Sandy argiles, humid and 
cultivated 

54 

A6 B6 C4 
57 59 57 

traditional urban area 

55 

A4 B4 C3 D4 
51 53 55 25. 

sandy argile humid area 

5 6 

A5 B5 C5 D5 
53 55 57 26 

sandy argiles 


It can be seen that the coding boundaries are adapted for the 
classification into subgroups of groups 2, 4, and 5, so that in each 
group the transformed channels will have a taxonomic importance for 
the binary image. 

8.1.2 Legend for the FOR I Map 

The legend for the printed FOR 1 map which is given in the appen- 
dix gives the main features and the dominating objects for them, with 
consideration of the spectral statistics of each classification group 
This legend is described in the following table. 


Table 8. A Legend of Map FOR-1 



Class 

Dominant 


Color 

Numb er 

Objects 

Features 


Spectral 

Signatures 

4 ['5 [ 6 | 7 


50 50 35 11 






shallow puddle 

puddle water s vege- 
tation of puddles, 
drled-up river tri- 
butaries 


puddle water, , , « 

vegetation ^ 10 


Miles i and Kousseri high chloro- 
forests, high chloro- phill vege- 
phyll vegetation tation 


43 39 49 24 


high chlorophyll 
vegetation, water 
marsh areas 

berber soil, culti- 
vation, vegetation 

modern residential 
urban center, vege- 
tation areas 


chlorophyll . „ . 

vegetation 42 35 45 21 


berber 


48 48 49 22 


green village 45 44 50 24 


humid argile soils, sandy argiles 49 49 53 25 

moderate chlorophyll and vegetation 


inhabited close plot traditional 58 62 60 27 

regions urban area 

dense plot inhabited traditional 57 59 57 26 

areas urban area 


suburban, bare ground bare sandy 60 65 64 29 

argiles 

high reflection on sandy argiles 57 60 6l 28 

humid ground 

river terraces sandy argiles 53 55 67 26 

terraces and creeks sandy argiles 51 53 55 25 


sandy riverbanks , 
creek boundaries 


sand^ 


66 73 71 31 


zones projecting 
from the nagas and 
river edges 


sandy argiles :g g 7 g6 
nagas J 
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The classification in steps, using successive cutting off of groups J 
(spectral coassistance) as a function of their taxonomic importance, is | 
based on the analysis of transformed canal boundaries and their inter- I 
pretation on the binary image. I 

With respect to FOR-1, which occurs at the end of the rainy season 
in November, 1972, the image treated by FOR-2 observes the area in the 
middle of the dry season in February, 1973. 

By comparing the two classifications, one can very clearly see a 
reduction in the rivers and the change in the flooded and humid areas 
and the chlorophyll areas into dry areas. 



channel . 


In this step of the processing of the data, on the binary image we 
select the level intensities corresponding to a large number of fea- 
turcc The study of the intensity levels allows a first identification 
of the features. 

The first result of automatic classification into three boundaries 
which separates classes of equal population in Channel 4 is given by 
the image shown below (Figure 8b). One obtains the selection of the 
water and the vegetation associated with the water. 
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Figure 8.B. The FOR-2 ima e of Channel 4, rivers, and puddles. 


/nc 

In the channel *1 (MSS 7), reflection leve from 1 to 19 are assc-~ — 
dated with the following objects within the ^ame feature: 

- rivers 

- ponds and puddles 

- flooded areas 

- airport runways 

- very humid and high chlorophyll zones. 

The average luminescence of objects on this feature varies between 
0.73 to 13.87 W*m sr . It can be seen that the airport runway has 
a low luminescence level. 

We studied each intensity level before proceeding with the first 
classification. We observed trees near the transformed canals, and 
the resemblance oetween the isodensity levels of the four channels. 
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8.2.2 First Classification of FOR 2 into Seven Groups 


. y_ 


8.2.2.J, Aggregation Control Parameters 

Groups [ Transformed channels | 
and reflection levels 


Objects 


Boundaries 


SS 4 31 37 41 43 46 48 52 56 60 113 


SS 5 25 37 42 46 49 52 55 60 63 6? 127 


22 34 42 49 -54 58 63 [72 127 


15 19 21 23 25 27 29 33 I 63 


Features 


Cl D 1 

19 5 

Chari and Logone 
rivers 

BIO A9 C9 BIO 
65 81 79 36 

sandy banks, sandy 
river edges 

B1 C3 B3 A 2 

Miles! and Kousseri 

40 39 18 40 

wood 

C8 D9 

bare urban areas. 

66 31 

grediya 

B9 A7 

Hwy 40, diguel. 

63 54 

Sabangali, Kousseri, 


village, Chagoua 

BO B7 

dense plot urban 

26 56 

area 

C4 

dense plot urban 

49 

area and central 


urban area ' . 


rivers 


bare and dry 


vegetation 


argiles 


suburban and 
urban areas, 
sandy argiles 


dry sandy 
argiles 


dry sandy 
argiles 


The binary image of the group 2 shown below is an example of the 
st FOR-2 classification. 
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Figure 8.C. Binary image of Group 2 , 
is clear* on the image. 

T.ViN-:s t r*} [504 ( G) Columns i 

po.?? 


identification of sandy areas 
S59 < 0 > Intensities o n ass 


Studied Class. 2 


8.2.3 The second FOR classification into fifteen classes. 

In the second step, Grups 3 5 6, and 7 are subdivided as a function 
of the selection based on the proximity scale of reflection levels 
which characterize the main objects of taxonomy. 


8 . 2 . 3 . 1 Analysis of the Binary Image of Subdivided Group 3 
a- A ggregation control parameters 


l sub 

: groups 

' transformed canals 
and reflection levels 

t * 

Objects 

Features 

\ 31 

A1 B1 C4 D5 
38 36 44 22 

Milesi forest, 
Kousseri forest, 
I.Z.V.A.C. gardens, 
green village, old 
meanders 

high chlorophyll 
vegetation 

t ^ 32 

A4 B4 C3 D3 
42 43 41 19 

cultivated areas, 
bottom land culti- 
vation, overgrown 
bottom land regions 

chlorophyll 

vegetation 
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33 


Cl 3)1 
30 11 


riverbank water , 
puddle wet zones 

center of puddles 


water and 
humidity 

puddle water 






i 



34 


C2 D2 A3 B3 
38 17 40 39 


On the following page, two illustration of the binary image of 
group 3 are given. 


b- Coding boundaries for group 3 


Channels 



Boundaries 


Number of 



1 

2 

3 

4 

5 

6 

boundaries 

MSS 

4 , 

31 

38 

40 

41 

53 


5 

MSS 

5 

25 

36 

39 

42 

53 


5 

MS 3 

6 

5 

36 

39 

42 

! 59 


5 

M3 3 

il 

7 

15 

17 

19 

22 

63 

6 



«• -v:s 3S4 a so * . i o) channel i a 259 c 

F0R3 

class studied 


0 > intensities 9 a ssp 
31 



Figure 8.D. Group 3 binary image: Class 31. Identification of the 
high chlorophyll vegetation areas is very clear. 


r„*. a bim c 01 channel 1 a £!•=? < <? ) ant eh si ties o a sss 

_ 

class studied 33 



Figure 8.E. Binary image of group 3, class 32. Identification of 

average chlorophyll areas in relation to the puddles and 
the bottom land areas. 



a- Aggregation control parameters 


Subgroups Transformed channels Objects ■ Features 

and reflection levels 
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b- Coding boundaries for Group 6 


Channels 


Boundaries 


Number of 
Boundaries 



1 . 

2 

3 

4 

5 

MSS 4 

40 

46 

49 

52 

100 

MSS 5 

25 

! 52 

56 

60 

127 

MSS 6 

! 2 ' 

53 

55 

57 

| 106 

MSS 7 

1 

24 

26 

28 

34 


The study of the group 6 subgroups shows the overall drying pat- 
tern of the soil and of the vegetation of the terraces and the urban 
zone. The chlorophyll effect seems to disappear. It is possible to 
distinguish group families of soils, and the urban morphology, depen- 
ding on the construction types. 


•. n?WL i 

li&XiZS 35-* A 50-4 < O) channel ± & 259 < a ) intensities e a sss 

. F0R2 

class studied S 2 


I n . .. ■ f 

- ■ -vv,. .** ' , f 

f * * ~ ’ ' r 

, j- . ' .* ^ - .f 

r ffv’ -’v /. ‘ . • l 


/ r i 
.4 .*• 






-I T 'I 1 • 1 I — 1 1 1 J 1 ' 

150 200 250 


Figure 8.F. • Binary Image of Group 6: Class 62. Identification of 

several urban islands with high reflection. 


8. 2. 3*3 Analysis of the binary image of subdivided group 7. 
As examples of this analysis, we give Figures 8H and 8G: 




a- Aggregation control parameters 



Subgroups Transformed Channels 

Objects 

Features 

and Reflection Levels 



Traditional 
Urban Areas 


Urban areas with Modern urban 

high reflection, area, bare 

modern construction, argiles 
dried sandy argiles 







■ • I r _ - 

in sc>4 c o) channel i a as*. < «• ) intensities c A 255 

FOP ' 5 

'class' studied ' ?3 


% * S* 

- t/ 


; ; f 

* r 


( 5 •C- - ’ - Lu _ 

j . «.r % 

I * - * ev 1 

- » : _ *• ^jt> _ 

_.{ /“ 

. t , . -M»* . . «“ Mj. a " *■*“• *" 

i_. ■ * * -* v '*’■ % . 

"? r » . 5 s’ _ * : ■» 

- v. *i v«jt p . •«. 

_ £i^w rvr: i .. .• - • 
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100 ISO 
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Figure 8.H. Binary Image of FOR-2 : Drying out of bottom lane, river 
beds 3 and green urban areas. 


Subgroups 

f Transformed Channels 

and Reflection Levels 

Objects 

Features 

/ 

73 

' A1 B1 

41 • 44 

Dried-out puddle 
areas, urban areas, 
during drying out 

Dried bottom 
lands, urban 
areas 

74 

Cl D1 

47 22 

Dried puddle edges 

dried bottom 
lands 


























£SS3ff$g~Bma»« 


A study of the sub-groups of group 7 shows the drying out of humid 
argile regions including the urban areas. The urban areas differ ac- 
cording to modern or traditional high-density construction methods. 
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8.2.4 FOR-2 Classifications 




f. 

i 

8.2.5 L egend for the FOR- 2 Map 783 


The final 1 OR- 2 map legend shown in the appendix can be used to 
locate the main features and the dominating objects with consideration 
of the spectral statistics of each classification group. 


Color 

Class | 
Number 

Dominating 
Co.} ects 

Features 

Spectral 

Signature 

Blue 

1 

Chari and Logone rivers: 

rivers 

40 

33 

19 

5 


33 

river edge waters, wet 
puddle areas 

water and 
humidity 

41 

39 

30 

11 


92 


Green 


Red 


Chestnut 

color 


Yellow 


34 

center of puddles 

puddle 

40 

39 

38 

31 

Milesi and Kousseri 
forests, IZVAC gar- 
dens, green village, 
old meanders 

high chloro- 
phyll vege- 
tation. 

38 

36 

44 

32 

cultivation, bottom 
land overgrown areas 

chlorophyll 

vegetation 

42 

43 

41 

72 

high reflection urban 
areas, modern con- 
struction, dried ar- 
gile soil 

modern 
urban area, 
bare argile s 

45 

49 

53 

71 

Traditional plot 
areas, sandy argile 
soils 

traditional 
urban areas 

47 

51 

50 

63 

traditional residence 
urban areas, dried 
sandy argile areas 

traditional 
urban areas, 
sandy argile 
areas 

50 

56 

54 

5 

urban and suburban 
traditional area, 

Hwy 40, Diguel, 
Sabangali, Kousseri- 
village Chagoua, ... 

Traditional 
urban area, 
sandy argile 
area 

54 

63 

61 

61 

levels of terraces, 
river edges 

dried sandy 
argiles 

50 

58 

58 

62 

urban islands with 
high reflection, 
terraced rivers 

dried sandy 
argiles 

46 

53 

56 

73 

dried edges of pud- 
dles, urban zones 
undergoing drying 

dried bot- 
tom lands , 
urban areas 

42 

44 

47 

74 

dried edges of pud- 
dles 

dried 

bottom lands 

45 

48 

47 

2 

bare and dry sandy 
banks revealed by 
settling 

sandbanks , 
sandy river 
edges 

65 

81 

79 

4 

sandy river edges 
and halomorphic 
soils 

urban bare 

areas 

nagas 

58 

68 

66 


17 


22 


19 


26 


23 


25 


28 


27 


27 


23 

22 


36 


31 




toil 


m 

s 
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8.3 FOR- 3 PROCESSING SEQUENCE 


/84 
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The classification in steps, using successive cutting off of 
groups (spectral coassistance) as a function of their taxonomic im- 
portance, is based on the analysis of transformed canal boundaries 
and their interpretation on the binary image. 

The following tables illustrate this for FOR-3 

8.3.1 First Classification Gf FOR-3 into Five Groups 

8. 3. 1.1 Control Parameters of Aggregation Before Co-assistance. 


Groups 

... , 

Transformed Canals 
and Reflection Levels 

Taxonomy f 

Post-Assistance 

1 

CIO D9 A9 B9 
62 28 53 64 

sand 

2 . 

Cl D1 

18 4 

rivers 

3 

B2 A2 D3 

34 36 16 

vegetation, puddle 
waters 

.'4 

1)8 A7 B7 
24 ‘ 47, 54 

sandy argiles 

5 

B3 A 6 C5 

46 42 45 

urban, argiles 


8. 3.1. 2 Separation coding boundaries 


Channel ' 

H 


3 

4 


■HI 

1 6 : 

9 

1 10 


Number of 
Boundaries 

KSS 

a 


33 

36 1 

38 j 

40 

42 

- 



52 | 

m 


10 

MSS 

5- 

Wk 

31 

36 

42 

45 

48 

1 

57 

62 

104 


10 

MSS 

6 ' 

14 

~ 

30 

37 

40 

| 45 

! 48 

.... 

51 

54 

57 

106 

11 

.tuO 

7 

2 

: 6 



18 

20 

! 22 

24 

26 

63 


10 





8. 3.1. 3 Division of this Classification: Example of Group 3 

In this first step, the groups 3, 4, and 5 are divided up, based 
on the selection of transformed channels as a function of their fea- 
ture significance. • 
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Therefore, group 3 is subdivided into five subgroups, 31 to 35, 
in which the selection of channels Bl-.and B2 and A2 characterizes: itT 
particular the subgroups 31, 32, and 33, where the high-chlorophyll 
vegetation dominates. The selection of channel D3 characterizes sub- 
groups 34. and 3^, where stagnant water of puddles dominates. 


The following tables illustrate this classification of Group 3. 
Sub Groups I ' Transformed Channels 1 Taxonomy 


' Transformed Channels 
and Reflection Levels 


31 

A1 B1 
32 29 

• 32 

C 6 D 6 


42 21 

33 

I>4 


15 • 

34 

Cl D1 


27 10 

35 

A 6 36 C5 


38 38 36 


Taxonomy 

Post-Assistance 

highly-wooded chlorophyll 
vegetation, humid ground 

High chlorophyll vegeta- 
tion, woods, and shrubs 

Puddle vegetation in humid 
areas 


Pixels of transition between 
water and riverbanks 


Low Chlorophyll humid zones, 
bottom lands. 


The coding boundaries were modified, as given in the following 
table : 


Channels 


Boundaries 


] 

1 

2 j 

3 

4 

5 

6 

7 

MSS 4 

28 

33 

34 

35 

36 

; 38 

1 

46 

MSS 5 

22 

“ 

30 

32 

33 

35 

37 

49 

~MSS 6 


~ 

32 

33 

35 

38 

49 

MSS 7 

5 

12 

14 

15 

17 

19 

26 


Numbers of 
Boundaries 


8. 3.1.4 First FOR 3 classification into 16 classes 
The groups 4 and 5 are subdivided using the descending method. 
Finally, one obtains a first classification into 16 classes (five 
groups), which is summarized in the table on the following page. This 
first classification will be completely discussed in the following. 




Groups 

Transformed 

Channels 

Taxonomy 

/ 

1 

CIO, D9 r A9, B9 

Bare sand, sandy 
argiles 

2 

Cl, D1 

Rivers 


31- 

Al, B1 

Forests, humid soil 


32 

C 6 , D 6 

Forests, shrub carpet 

3 

33 

D4 

Humid area and puddle shrub 
vegetation 


34 

Cl, D1 

Transition pixels 


35 

A 6 , B 6 , C5 

Humid low chlorophyll bottom 
land areas, puddle waters. 


41 

A4, B5 

Naga soil 


42 

. D4 

Naga soil 

4 

43 

C3, B4, B3 • 

Sandy Argiles and urban 
areas with high density 


44 

Cl 

Sandy argiles and urban 
high density area 




Modern urban areas and 
vegetation drying out 


51 

Al, B1 

5 

52 

Cl, D1 

modern urban areas and 
vegetation drying out 


53 

A2 , B2 

very high density urban 
plot area 


54 

B 6 , C5, D5, A5 

low terraces 


55 

A4, C 2 

high reflection urban area 


8-3.2 Second classification of FOR - 3 into nine g roups. 

Based on successive classification into five groups discussed 
above, a second classification is then carried out which leads to the 
final legend for Map FOR-3. ... , 

This second classification is a result of the first one with the 
following intermediate step: 































































Class 


Transformed Channels 
and Reflection Levels 



CIO 

D9 

A9 

B9 

62 

28 

54 

65 

Cl 

D1 


■ .*<: " 

17 

4 



B1 

B2 

A2 



34 

36 

* 

D3 




13 




B8. 

D8 



58 

26 





Obj ects 


bars and river 
edges 


rivers 


forests, bottom land 
vegetation, green 
villages 


puddles, dry river 
beds 


nagas 


African urban area 
with low density 
plots 


green humid berber 
areas, partially 
'urban (no large 

trees*)' 


river edges and 
terraces, low den- 
sity urban areas 


modern urban central 
area, traditional cen- 
tral urban area with 
high density 


Features 


sand 


water 


vegetation 


water 


argiles 


argil es 


humid 

argiles 


argile 
silts . 
a 

argiles 


urban area 


Aggregation control .parameters. 

The analysis of the binary image of the 9 groups allows one to 
introduce a subdivision of Groups 3 a 8, and 9. Identificationof the 
dominating objects of the features is done from the intermediate out- 
puts of the binary image during the spectral coassistance step. 

The tendencies of the spectral statistics of the groups confirm 
the identification of the features. 

























































Transformed Channels 
and Reflection Levels 


Oojects 


Features 



" C3. D4" Al 331 ' 
41* 21 33 30 

Milesi and Kousseri 
forests, IZVAC, cen- 
tral green modern ur- 
ban area, Chari dried 
tributarites , green 
zones of cultivated 
puddles 

C2 D3 A2 B2 
.37 18 34 32 

vegetation on argiles 
of bottom lands and 
puddles . 

1)3 A/j R4 t . 

11 40 38 

small islands, un- 
covered, shrub vege- 
tation on sand-silt 
background. 


vegetation 
forests, dense 
urban areas. 


33 16 36 34 


shrubbery on bottom 
lands, acacia supal 


High chloro 
phyll vege- 
tation . 


vegetation of 
bottom lands 


bottom land 
vegetatidn. 


Separation coding boundaries. 

The selection of the transformed channel boundaries is done as a 
function of the feature significance. 

Parameters for controlling aggregation. 


Groups 


Channels 



39,98 

1,75 


-29,99 

2,64' 

40,24 

4,94 

32,27 

2,90 

37,09 

1,63 


im 


Spectral statistics of groups 


The tendencies of the spectral statistics of the subgroups confirm 
the feature identification. 

















































Spectral Statistics of Groups. 

The method utilized is the same for the division of Group 3 
8. 3*2.4 The 9 Groups: Subdivision into two subgroups. 
Channels 


SS 6 


Boundaries 
~ 3 


30 i 

39 

22 

39 

14 

42 ! 


43 

60 

45 

48 

48 

108 ] 


Number of 
Boundaries 


1 . A2 . 

B2 


! ; 39 V 

42 

. • i 

C2 

D2 

r 

43- 

20 ; 



Modern central urban urban area 
area 3 dry argiles of 
bottom lands 3 dry ar- 
giles of low terraces, 
vegetation in sandy 
soils 


traditional central 
urban area with high 
density 


urban area 


Parameters for controlling aggregation 


Grouos 


Channels 


1 

9 

MSS 5 ' 

MSS* 6 

MSS 7 

91 39,24 

41,72 

45,29 

21,77' 

' 1,25 

2,63 

2,07 

1,29 

92 41,68 

45,09 

43,44 

19,76 

1,54 

2,13 

1,69 

1,02 
























Spectral Statistics of Groups 


The method used is the same as for the division of group 3. 


8.3.3 Legend for Map FOR- 3 || 

The map legend FOR-3 3 which was printed and presented in the appen- 
dix, allows one to locate the various features and their dominating ob- 
jects, with consideration of the spectral statistics of each group of 
classification. 


Color 


Blue 


Number of 
Class 



Green 



Dominant 

Objects 


Chari and Logone 
rivers, low tur- 
bidity water 

puddle water, 
flooded areas, 
dried up tri- 
butaries 

Milesi and Kou- 
seri forest, 
green village, 
shrub cover 

Dense shrub of 
flooded and pud- 
dle areas 

bottom land vege- 
tation and shal- 
low water 

bottom land 
vegetation and 
low terraces 


Features 


Spectral 

Signature 


rivers 


39 31 


puddle water 34 32 


high chloro- 33 30 

phyll vege- 
tation 

Chlorophyll 34 32 37 

vegetation 

vegetation and 

water j 40 38 30 


vegetation 


modern urban area humid and 
and humid chloro- urban areas 
phyll areas 


modern urban 
areas and river 
terraces 


urban area 


traditional urban I urban area 
area of high den- 


34 34 


39 39 


42 45 


42 45 43 


argiles 






Table 8.J. Legend for Map FOR- 3 
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8.4 Sequence of FOR-1.2 Processing, or FOR 1 and FOR 2 Diachrony 

We decided to perform st . diachronic treatment by superimposing the 
five measurement channels 3 that is, the four channels MSS- 4 , MSS-5, 

MSS- 6 3 MSS-T of the data from November 1, 1972, (FOR-1) and the four 
channels MSS-43 MSS-5, MSS-6 3 MSS- 7 'from the data from February 17, 

1973 (FOR- 2) . 

8.4.1 The Classification Method in Steps 

The classification method in steps assumes that one can perform 
a spectral character assistance procedure and a feature procedure. 

In the case of a Aiachronip<; study 9 the spectral statistics of the 
groups are particularly Interesting factors to be taken into account 
during the interpretation. From the difference of the forms of the 
spectra of the two days 3 one can appreciate the stability of the feature 
importance of one group over time. This means that one can detect changes 
in the landscape over time. In the case of LANCHAD, the changes are 
essentially determined for the variation of the humidity of the ground 
and the vegetation functions over the year. 

We will give three examples of this type of interpretation in Figure 


I. 


§ 

3 

Ft 





l 

I 

•A 

% 

i 


i 

I 


-i$ 

* 

4 

t 

% 


•I 


104 





8.K. In the first case,, we are dealing with dense forests, (Group 31), 
whose spectra essentially have the same form for either observation 
date. The chlorophyll characteristic is dominant and is relatively con- 
stant over time. In the second case, we are considering chlorophyll 
background (Group 33) whose vegetation function changes from November 
to February. One can see that the [j^fiectidtespectrum of February is 
slightly flatter than in November, which translates into a certain drying 
of the vegetation. In the third case. Group 51, the variation between 
the two dates is even more distinct. We are dealing with bare and dry 
sand banks in February, but they are covered with water in November at 
the beginning of the settling of the Chari and Logone rivers. 


The systematic observation of the spectral signatures, therefore, 
allows us to take into account a diachronic study in the interpretation 
of the kinematic characteristics. 




• November./ 1972 


f^.f;February 1973 


mm mm mm Has? hs m hsss 

a- Spectral stability of dense forests (Group 31). 


HSS7 



November'; 1972 


i.ary, 1973 


BSS5 HS3S HSS? NSS5 KSSS 1 H5S7 


b~ Drying of chlorophyll bottom lands (Group 33). 






TO5 


c- Rough spectral variation produced by the drying of certain sand' 
banks (Group 51). 

Figure 8.K: Three examples of spectral statistics of the FOR 1.2 

v diachrony . ' ■ 
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8.^.2 Classifications of FOR-1.2 


First Classification: 


' /91 

Second Classification 


35 groups and 
subgroups 


First Map C 3/77) 
23 groups and 
subgroups 


Final Map (8/77) 
31 groups and 
sub groups 


•21 + 25 

21 

2Z’^:.y v 

22 

.23 i>24 • ' f - 

23 + 24 

* ' > " " :'V 

f 

25 

31 

31 

.32 + 35 

32 + 35 

*33 i 34 

33 + 34 

• c ■ • . , " 












8.4.3 Legend of Map FOR 1.2 


• ^ Tfae final legend for Map FOR. 1.2 given in the appendix allows one 
to -identify the,- various features and tile ir main objects, with consider- 
ation of' the special statistics of each- group- of classifications. This 
legend is described in the following table: 

Table feh':. Legend of Map FOR 1.2 ■■■■■■■-*• : 


Spectral 
di achronic 
jgodr: -channel 
s’ignature 


Color 

Class 

Dominant Obj ects 

Features 


Numb er 

V- 5V ; . . 

/• 


Blue 


Green 


42 deep rivers, low tur- 

bidity 3 on the bottom 

45+46 shallow rivers decrea- 
sing turbidity toward 
the bottom 


47 -river banks and puddle 
i centers during .drying 

44 + 48 puddles drying up 

4l vegetation on puddle 

waters during drying 

43 argile depression, 

vertiso ls h um id t o 

hi gh 

t ent ■ /Vf v ' t;v/ " v ' •••• 

31 Milesi woods, Kousseri 

woods , IZVAC gardens , 
parks, forests 

32+35 high chlorophyll puddle' 
area shrubs, ponds 

33+34 bottom lands and old 

tributaries, very humid 
and high chlorophyll, 
cultivated 

1 flooded areas drying 

up with dry vegetation 


puddle 

water 


vegetation 

argile 1 

bottom 

lands 

high chlo- 
rophyll 
vegetation 


phyll ve- 
getation 

chlorophyll 
bottom In- 
lands 

depres- 
sions and 5- 
flooded 
zones 


31 

11 

36 

12 

44 

12 

37 

17 

34 

19 

39 

23 

36 

25 

39 

22 

43 

20 

60 

21 


Red 


Chestnut 

color 


Ik 

Humid argile and chloro- 
phyll areas drying up 

terraces 

43 

22 

Residential areas, green 

urban 

45 


villages, drying tendency 

area, green 

23+24 


village 


flooded areas to planted 
areas 

flooded 

areas 

52 


72 

urban areas and islands 

modern 

55 


with modern construction 

! urban 




areas 


71 

traditional central high 

traditi- 

60 


density area 

onal urban 



area 


82+83 

traditional low density 

suburban 

64 


areas and suburban areas 

area 


92 

suburban low density 

suburban 

62 


area with bare areas 

area, 

sandy 

argiles 


61 

river banks and terraces 
to cultivated areas 

river bank 

s 

66 

62 

riverbanks and creeks 

nagas and 



and nagas, to cultivated 

riverb anks 

72 


areas . 



63 

urban areas and river- 

bare ar- 



banks to dry areas, sub- 

giles. 

64 


urban areas 

sandy areas 

91 

terraces and riverbanks 

sandy 

argiles 

49 

93 

riverbanks to dry areas 

sandy 

argiles 

56 

73 

suburban and riverbanks 

sandy 

argiles 

55 

75+76 

riverbanks and cultivated 

sandy 

53 


terraces 

argiles 


21 

riverbanks and cultivated 

terraces 

49 


terraces 

and river 

- 



banks 


25 

argile sandy areas to 
dry areas 

terraces 

49 




s.'S 


m 


H 

I 


ts 


Yellow 


81+84 


51 


52 

53 

54 


bare and dry areas 3 sub- 
urban areas 5 (Miles! 3 
Chagoua 3 Denib e) 

Dry and bare sandbanks 
on the second date 3 
and covered by water 
on the first date 

banks and riverbanks 
uncovered by settling 


very dry and bare sands 
on both dates 

sandy riverbanks and 
halomorphic soils 


suburban 
areas., be 
spaces 

banks and 
sand 


sandy ar- 

gile 

silts 

sandy 

argiles 

nagas and 
riverbanks 


_ J.O wicumei i^; 3 ^-5 on the second date is 

channel (6). MSS-7 of the first date is channel (4), and MSS-7 on the 
second date is channel (8). 
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~ andloR-i^ ° f LSKendS ° f Classi f ^ations FOR 1, FQR-2, P0R 3 

This correspondence will be analyzed in a subsequent study, which 
will take into account the truth-terrain data. 
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CHAPTER 9 - GEOGRAPHIC AND CARTOGRAPHIC INTERPRETATION 

The geographic interpretation of the results of classification for 
the maps presented here should be essentially based on the specific 
features of the landscapes observed on the recording dates by the 
SAT satellite. , . / /.. . . 

The detailed geographic analysis of the results trill be the topic 
of a separate publication, and there will be a Certain delay with res- 
pect to the numerical analysis. This will allow a comparison of the 
detailed results, and the totality of the information collected during 
ground truth surveys. 

We would like to stress that for the Sahara-Sudan j^n^^^nment^^ 
the variation of the surface phenomena related to the water cycle and 
to the vegetation cycle is extremely rapid. During the two observation 
dates by the LANDSAT satellite, for example, between November 1, 1972, 
and February 17, 1973, one finds completely opposite times during the 
water and vegetation cycles observed on the ground. Considering the 
orbital variations, the pixel density (80 m x 80 m) will probably be 
difficult to superimpose between the two dates. The variation of the 
spectral signatures of the pixels is a complex phenomenon related to the 
size of the observed objects and the increasing or decreasing diffusion 
of objects located within the pixels. 

A few specific features can be described as an example;... 


9.1 Traditional Urban Area Feature 

The traditional poto-poto habitation feature can or cannot in- 
clude the following: 

- tile roofs; ’ 

- a tree in the yard 

- a larger or smaller yard 

- a poto-poto enclosure wall 

- cases with seco (thatch) roof. 

- !■ : 




The traditional habitation islands can be limited by very narrow 
streets or can include areas of habitation which are widely dispersed 
over bare ground without enclosures or fences. This is what we will 
call traditional habitation areas with high density or with low density 


9.2 Modern Urban Area Features 

Modern constructions for habitation or other purposes can be 
constructed in high-density areas or in parks and gardens. 

The streets can be tarred or can be simple tracks. The avenue 
can be bare or have tree borders . The tiled roofs or concrete terraces 
can be partially covered by the branches of large trees, or can reflect 
greatly from bare ground. 

In the two cases, traditional areas or modern areas, the soil can 
be saturated with water or can be completely dry. 
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9.3 Rural Landscape Feature 


These features are known from pedological maps, or vegetation maps, 
witif scales between 1:200:000° to 1:500,000°. In these maps, division into 
zones results in two generalizations . The aerial surveys available come 
from different dates. 


The results of the classification of the numerical LANDSAT satellite 
data a very detailed approach to this type of environment; for 

example, large variations in the ground reflection were observed, depen- 
ding on the variation of the cultural cycle. The chlorophyll effect 
can dominate, depending on whether the water layer in the ground has 
already (or not yet) allowed the appearance of shrub vegetation. It also 
depends on the seasonal cycle variations of natural vegetation, and 
whether or not this leads to a detectable chlorophyll effect. The 
variation of the chlorophyll effect is very rapid for the following, 
for the Sahara-Sudan environment: 


i 


the forest galleries of the berger 'riverbanks, 

the shrub and pine foliage 

low density acacia or mimosa areas, 

low areas covered with shrubs with t;?^(to|iogpriepupv vegetation 
burned-over cultivated areas 

r ‘ 

rainy cultivated areas 
pond cultivated areas. 


at the scale of Territory Management 


The cartography presented here on the scale of 1: 50, 00ti° and 
1:100,000° shows the possibility of working at scales which are useful 
for managing territories, using numerical data from the LANDSAT satel- 
lite. This thematic analytic cartography must be followed by thematic 
cartography which is more synthetic and represents the step of geographic 
interpretation of results. This is done at a separate time. 

/IQ- 

We can only present the results of this synthesis work when the 
analytic cargotraphic works for the three 1975 and 1976 images appear for 
the same zone (FOR 4, FOR 5, FOR 6). 


The synthesis will therefore be the result of multiple season and 
multiple year analyses which we have performed. We believe that this 
is the only way of performing research in this area. 




CHAPTER 10 - LANCHAD PERSPECTIVES FOR THE PLANNED ANALYSIS OF THE 

SAHARA-SUDAN ENVIRONMENT 
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The LANCHAD atlas is one of the important tools for the analysis 
of the natural geographic environment or the cultivated environment 
in an area which extends from 7°N to 24°N and 12°E to 25°E. This is 
approximately 2 , 500 , 000 km 2 . 

This project has been presently proposed to the assistance and 
cooperation fund by the Tchad ministry of the economy. It is planned 
that a reference document will be prepared for managing the territory. 

The types of terrain which are planned to be analyzed are the fol- 
lowing, within the framework defined by these six parts of the regional 
Tchad atlas: 


Part 1: Bardai Oasis, Pic Tousside, Emi Koussi; Faya Oasis, 

Ounianga Oasis, Sahara region. 

Part 2: Mao Oasis and Diwe basins, Ouadi cultures, pastures; Liwa 
and sodium depressions; 

Bol Polders and Chari delta; 

Sahelia region. 

Part 3- Djamena zone, urbanisation and farms; deforestation due 

to urban economy; Bongor and rice plantations; Sahelia-Sudan 
area 

Part 4: Lai and Tandjile rice plantation, Sategui-Deressia; Moundou, 
Gounou Gaya, cotton fields; Pala,. Fianga, Lere, mining area; 
Goundi, Dono Manga, Dob^, Koumra, les Koros, le Mandoul; 
the tropical Sudan area'. 

Part 5:- Sarh, Moissala; Am-Timam, Haraze Mangueigne, le Salamat, 

cotton plantation, sugar cane, peanuts, explotetion of the 
the Karite region 

Part 6: Abeche, Ati, Mongo; peanut plantations, ouadi cultivation, 
pasture. 


In each part, the dynamic analysis of the area will be. directed by 
the condition of providing the useful and new information, capable of 
supporting maj.or management decisions of the territory within the gene- 
ral project of development of the Tchad up to 1980. 






The following topics will be discussed: 
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- urbanization and effect of urban economy on the area, 

- traditional and modern agriculture, 

- agriculture and industrialization, cotton, rice, peanuts 

-- arborculture and forestry 

- erosions of the ground 

- geology and geomorphology 

- hydrogeology and hydrology, 

- study of regions and territories 

- use of various soil areas. 

Considering the remote detection methods used, it is predictable 
that their number and their capacity will overlap. 

As far as the wavelength is concerned, we would like to state that 
in 1978 the LANDSAT-D satellite will allow operation in the thermal 
infrared range (10 to 12 M) with a resolution comparable to that of 
the present LANDSAT satellites. This new possibility will allow a 
detailed study of the evaporative transferation of the vegetable growth. 
This is an essential factor for a good definition of cultivation. 

As far as the repetitiveness of the observations in time is con- 
cerned, there will be a substantial change in 1978 because of the 
launching of the European METEOSAT satellite. This will allow double 
observations in the visible and infrared ranges, with a resolution of 
about 1 km with a potential repetitiveness of 1 observation every 30 
minutes, day and night. This will mean that for the first time the 
Tchad will be observed during various seasons, in particularly, the 
rainy season. The accumulation of images during one day will allow 
the elimination of the massive cloud cover. The clouds are displaced 
rapidly enough during one day so that they do not cover the same part 
of the country over a period of 24 hours. There will also be substan- 
tial improvement in the areas of geology and hydrogeology, in parti- 
cular, the temperature differences on the ground as a function of the 
time of day and night will be valuable. The variation of these very 
differences will depend on a parameter called the thermal inertia, which 
is directly related to the characteristics of the rocks and the humidity 





The multiplicity of these new observations will mean that we wiill 4 

have to establish a reference data bank showing ground truth values f 

for each of the topics of this project. In this case, recourse to | 

the information is a requirement considering the complicated processing | 

of the data from remote sensing. 

LIST OF ABBREVIATIONS 

AEF - French Equitorial Africa 
APT - Automatic Picture Transmission 
ASECNA - Association for Airborne Navigation Safety 
CEMS - Center of Space Meteorological Studies 
CEPlCG - Study and Research Center for Geographic Cartography 
CIRCE - Interregional Center of Electronic Calculations 
CNRS - National Center for Scientific Research 
CTAMN - Remote Detection Center and Analysis of Natural Media 
DMG - Directorate .of Mines and Geology of the AEF 
EHESS = Advanced Social Science School 
ENSJF - Advanced Normal School for Young Girls 
ENSMP - National High School of Mines in Paris 
ENTP - National Public Works School 
EPHE - Practical Advanced School 
ERTS - Earth Resources Technology Satellite 
ESCG - Advanced Geographic Cartography School 
FAC - Aid and Cooperation Fund 
IGNF - National Geographic Institute of France 
IRC - Infra-Red Color 

MICA - Ministry of Industry, of Commerce and Artistry 
NOAA - National and Oceanic Atmospheric Administration 
ORSTOM - Scientific Research and Technical Office for Overseas 
Pixel - Picture Element 
RCP - Operative Research on the Program 

STISI - Data Processing and Statistical Industrial Service 
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